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Abstract 

We consider a system of non-relativistic spinless quantum parti- 
cles ( "electrons" ) interacting with a quantized scalar Bose field (whose 
excitations we call "photons"). We examine the case when the veloc- 
ity V of the electrons is small with respect to the one of the pho- 
tons, denoted hy c {v/c = e <^ 1). We show that dressed particle 
states exist (particles surrounded by "virtual photons"), which, up to 
terms of order (v/c)^, follow Hamiltonian dynamics. The effective A^- 
particle Hamiltonian contains the kinetic energies of the particles and 
Coulomb- like pair potentials at order (v/c)^ and the velocity depen- 
dent Darwin interaction and a mass renormalization at order {y/cf'. 
Beyond that order the effective dynamics are expected to be dissipa- 
tive. 

The main mathematical tool we use is adiabatic perturbation the- 
ory. However, in the present case there is no eigenvalue which is 
separated by a gap from the rest of the spectrum, but its role is taken 
by the bottom of the absolutely continuous spectrum, which is not an 
eigenvalue. Nevertheless we construct approximate dressed electrons 
subspaces, which are adiabatically invariant for the dynamics up to or- 
der {v I c) Y^ln[(z;/c)~-'-]. We also give an explicit expression for the non 
adiabatic transitions corresponding to emission of free photons. For 
the radiated energy we obtain the quantum analogue of the Larmor 
formula of classical electrodynamics. 

1 Introduction 

In a system of classical charges the interactions are mediated through the 
electromagnetic field. In the case when the velocities of the particles are 
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small with respect to the speed of light, it is possible, loosely speaking, to 
expand their dynamics in powers of v/c. 

The qualitative picture which emerges has three main features. Up to 
terms of order (f /c)^, the dynamics of the particles are still of Hamiltonian 
form. Moreover, 

• at leading order, (f/c)°, the retardation effects can be neglected and 
the particles interact through an instantaneous pair potential; 

• at order (f/c)^, the particles acquire an effective mass, due to the con- 
tribution of the electromagnetic mass and, to take into account the re- 
tardation effects, one has to add to the potential a velocity-dependent 
term (the so-called Darwin term). 

Including the terms of order (f/c)^, the dynamics are not Hamiltonian 
anymore, instead 

• there is dissipation of energy through radiation. In the dipole approx- 
imation, the rate of emitted energy is proportional to the acceleration 
of a particle squared. 

A formal derivation of this picture, which does not consider the problem 
of mass renormalization, can be found in |LaLij . A mathematical analysis in 
the framework of the Abraham model, i. e., for charges which have a rigid 
charge distribution, is given in [KuSpi] |KuSp2 . 

The above description is expected to remain true also for nonrelativistic 
quantum electrodynamics, where, neglecting the possibility of pair creation, 
one considers a system of nonrelativistic particles interacting with the 
quantized Maxwell field. In physical terms, the particles carry now a cloud 
of virtual photons, which makes them heavier, and interact exchanging them 
or dissipate energy through photons travelling freely to infinity. However, 
in quantum mechanics one describes the interaction of charged particles in 
good approximation by introducing instantaneous pair potentials and with- 
out treating the field as dynamical variable. If the particles move sufficiently 
slowly this is known to be a very good approximation. One goal of our pa- 
per is the mathematical derivation of quantum mechanics from a model of 
particles that are coupled to a quantized field, but do not interact directly. 

Instead of nonrelativistic QED we consider the massless Nelson model, 
which describes A^ spinless particles (which will be called "electrons") cou- 
pled to a scalar Bose field of zero mass (whose excitations will be called 
"photons"). In spite of the simplifications introduced, this model is expected 
to retain the main physical features of the original one. Therefore, since 
its introduction by Nelson |Nej . who analyzed its ultraviolet behavior, it has 
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been extensively studied to get information about the spectral and scattering 
features of QED, mostly concerning its infrared behavior (not pretending to 
be exhaustive, some papers related to this aspect are |Ar] [F?] |LMSj [Pi] and 
references therein). The recent monograph by Spohn |Sp| contains detailed 
descriptions of the classical and the quantum mechanical models and results. 

In this paper we define and analyze the dynamics of dressed electron 
states in the Nelson model in the limit of small particle velocities. Loosely 
speaking a dressed electron is a bare electron dragging with it a cloud of 
"virtual" photons. We show that the dynamics of dressed electrons have the 
features discussed above: an instantaneous pair interaction at leading order 
(v/c)^ and a renormalized mass together with the velocity-dependent Darwin 
interaction at order {v/cf'. We also provide an analogue to the classical 
Larmor formula for the radiated energy, i.e. for the energy carried away by 
"real" photons travelling freely to infinity (the heuristic concept of "real" 
and "virtual" photons will be made precise below). It is important to note 
that we consider the massless Nelson model with an ultraviolet cutoff but 
no infrared cutoff. Indeed, the leading order effective dynamics were already 
analyzed in |Te2| assuming an infrared regularization. One expects, and we 
will show it in this paper, that the dynamics of the dressed electrons even at 
higher orders are essentially independent of an infrared regularization, but 
the radiation is instead very sensitive to it, which makes the mathematical 
analysis much more delicate. 

To explain in more detail the kind of scaling we are interested in it is 
convenient to look first at the classical case. The classical equations of motion 
for particles with positions q^, mass ntj and a rigid "charge" distribution 
Qj coupled to a scalar field 0(a;, t) with propagation speed c are 

1 ^ 

mjq^it) = - J dx{V^(P){x,t)g,{x-q^{t)), 1<J<N. (2) 

We assume for simplicity that 

gj{x) = ej(p{x), 

where the form factor cp gives rise to a sharp ultraviolet cutoff, 

otherwise . 



^We use bold italic font, x, to denote vectors in M'^. The only exception to this, since 
there is no possibility of misunderstanding, is the three-dimensional momentum of the 
photons, denoted by k. The lightface font, x, will be used to denote vectors in M^^. 
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Taking formally the limit c — > oo in ([T]), we get the Poisson equation for 
the field, and so, eliminating the field from ([2]), we obtain equations of mo- 
tion describing particles interacting through smeared Coulomb potentials. 
Mass renormalization does not appear at the leading order. 

Instead of taking c — oo, we regard as more natural to explore the regime 
of particle properties which gives rise to effective equations. Indeed, if we 
look at heavy particles for long times, i. e., if we substitute t' = et and 
m'j = e^TJij in ([T]) and ([2]), we find that, up to rescaling, the limit e — > is 
equivalent to the limit c oo. After quantization however, the two limiting 
procedures are not equivalent anymore. 

The case c ^ oo was analyzed by Davies |Daj and by Hiroshima [Hij , who 
at the same time removed the ultraviolet cutoff, applying methods from the 
theory of the weak coupling limit. In this paper we analyze the limit £ — > 0. 

We define now briefly the massless Nelson model (whose presentation will 
be completed in section [2]), we state our main results and the principal ideas 
of the proof and then compare them to the above mentioned weak coupling 
limit. 

The model is obtained through canonical quantization of the classical 
system described by ([T]) and ([2]). The state space for spinless particles is 
L'^{E?^) and the Hamiltonian, assuming for simplicity that all the particles 
have equal mass, is given by (we switch to natural units, fixing h = 1 and 
c=l) 

1 ^ 

As explained above, we consider the case of heavy particles, i. e., 

m = 5"^ < £ < 1, (4) 
therefore the Hamiltonian becomes 

H; = -jA^,=:^p\ xeR'"", p:=-ieV.. (5) 

The particles are coupled to a scalar field, whose states are elements of 
the bosonic Fock space over L^(M'^), defined by 

^■■=®M=,®fs)L\R\ (6) 

where denotes the M-times symmetric tensor product and ®°^^L^(M^) := 
C. We denote by Qm the projector on the M-particles subspace of =^ and 
by VLy the vector (1, 0, . . .) G called the Fock vacuum. 
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The Hamiltonian for the free bosonic field is 

Hr.= dr{\k\), (7) 

where k is the momentum of the photons (the reader who is not familiar with 
the notation can find more details in section [2]). 

The particle j is linearly coupled to the field through the interaction 
Hamiltonian 

Hi,j := dy(f){y)gj{y-Xj), (8) 

where (p is the field operator in position representation. 

The state space of the combined system particles + field is 

^ := L^{R^^) ® ^ ~ L^{R^^, ^) (9) 

and its time evolution is generated by the Hamiltonian 

N 

:= H'q + J^Hij, i/o' -^p® l + l®dr(|fc|), (10) 

with domain 

J^o-= H\R^^,^)nL^{W^^,D{Hi)), H\W^^ , ^) := Dip" (g) 1) , (11) 

which is a Hilbert space with the graph norm associated to Hq. 

Note that there are no direct forces acting between the particles, all the 
interactions are mediated through the field. 

Our goal is to understand the dynamics of the particles for times of order 
e^^. It is necessary to look at long times in order to see nontrivial dynamics 
of the particles, because, since their mass is of order 0{6~'^) and we con- 
sider states of finite kinetic energy, their velocity is of order 0{e). However, 
since the coupling between the electrons and the field is not small, standard 
perturbation theory is of no use initially. Indeed, since the charge of the 
particles is of order one, the local deformation of the field, i. e. the "cloud of 
virtual photons" , is of order one. However, the influence of real photons with 
finite energy and momentum on the dynamics of the heavy electrons, whose 
mass is of order is small. Hence one expects that the coupling between 
properly defined dressed electron states and real photons is small. To make 
this precise, we construct a dressing transformation : Jif, which 

allows us to define the dressed particle states as follows. In the new repre- 
sentation the vacuum sector L^(M^^) ® Qp of ^ corresponds to states of 
dressed electrons without real photons, while in the original Hilbert space a 
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state of dressed electrons with M real photons would be a linear combination 
of states of the form 

(V^ ® a(/i)* ■ • • aifMT^F) with G L2(R3^) , /i, . . . , /m G ^^(MS) . 

Recall that Qm denotes the projector on the M-particles subspace of Fock 
space, then the projector on the subspace corresponding to dressed electrons 
with M real photons is 

In a nutshell, the main results we prove are the following: the subspaces 
are approximately invariant under the dynamics generated by for 
times of order e~^. For states inside such a subspace the dynamics of the 
particles can be described by an effective Hamiltonian for the particles alone 
on the above time scale and with errors of order log(£~^)) . Finally we 
can compute the leading order part of the state which leaves the subspace 
PqM' under the time evolution, which corresponds to emission of real pho- 
tons. The formula for the energy of the real photons traveling to infinity, i.e. 
the radiated energy, yields a quantum mechanical analogue of the classical 
Larmor formula for the radiation of accelerated charges. 

Before we can state our results in detail, we need to explain the adia- 
batic structure of the problem in some detail. The Hamiltonian is the 
perturbation of a fibered Hamiltonian, because, since Hi^ depends only on 
the configuration Xj of the jth particle, the operator 

N 

H^^{x):=dV{\k\) + Y,Hi,{x,) (12) 
acts on ^ for every fixed x G M'^^. This means that 

r® 

H' = ~-A,®1+ / fffib(x). 

Note the structural similarity with the Born-Oppenheimer approxima- 
tion. There the Hamiltonian describes the interaction between nuclei and 
electrons in a molecule and the former have a mass of order with 
respect to the latter (the typical spectrum of iffib(x) for a diatomic molecule 
is shown in figure [T]). In the present case, the particles take the role of the 
nuclei, and the bosons the one of the electrons. 

In contrast to the molecular case however, in the Nelson model if fib (a;) has 
typically (see lemma [2] and corollaries [1] and [2]) purely absolutely continuous 
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Figure 1: Schematic spectrum of the fibcrcd Hamiltonian in the case of a diatomic 
molecule for energies below the dissociation threshold (r = |£Ci — a;2|). The different 
eigenvalues are pointwise separated by a gap. 



ta(H.,(x)) 




Figure 2: Spectrum of iJfib(a;) for TV = 2 (r = |a;i — a;2|)- The spectrum is absolutely 
continuous, there is no eigenvalue at the bottom. The oscillations in r caused by the sharp 
ultraviolet cutoff are irrelevant, and therefore we do not show them. 
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spectrum, which does not display a structure with pointwise separated bands 
{absence of both eigenvalues and spectral gap) (see figure 

The bottom of the spectrum, E{x), can be exphcitly calculated, 

1 ^ 

Eix) = - ^l^.,(a^. - Xj) + Co, (13) 



where 



and 



Vij{z) := - dvdw — — — , (14) 

' .B3 4:n\v -w\ 



eo:=--yf (15) 



2p^iM3xR3 47r\v-w\ 

The effective pair-potential Vij (z) coincides, up to the sign, with the electro- 
static interaction energy of the charge distributions Qi and qj at distance z, 
while Co is the sum of all the self-energies. 

E{x) becomes an eigenvalue of Hfa,{x) if the total charge of the particles 
system is equal to zero, as it happens for example in presence of an infrared 
cutoff. In this case, it is possible to build for every x a unitary operator, 
V{x), called the dressing operator, which diagonalizes H{a^{x) in the sense 
that 

i^fib(x) = Vix)HiVixy + E{x) . (16) 

Exploiting this remark, we approximate the time evolution generated by 
if^ in two steps. First we define an infrared regular Hamiltonian, H'^''^, where 
the form factor if (equation ([H])) is replaced by 

^^(^)^^|(2^)-'/' cr<\k\<A, ^^^^ 
1 otherwise . 

Proposition, (see propositionUl) Let C{J^q^'^,J^) be the space of bounded 
linear operators from J^q^"^ to M' equipped with the operator norm. It holds 
then that 



e 



itHVe ^-itH^-'^/e\\ ^ r\ + \- (^ 



e 



where M'q^'^ := D[{HqY^'^^ with the corresponding graph 



norm. 



Choosing a as a power of e we can then replace the original dynamics 
with infrared regular ones. The latter contain however two parameters, e 
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and (7, therefore it is necessary to control the behavior of all the quantities 
that appear with respect to both. 

The advantage is that for H'^''^ we can build an approximate dressing op- 
erator '^e,ai acting on the whole Hilbert space which is defined for every 
positive 0", but whose limit when a — > 0^ does not exist if the system has a 
total charge different from zero. '^e,a is unitary and can be expanded in a 
series of powers of e, with coefficients which are at most logarithmically di- 
vergent in a. Moreover, the zero order coefficient is given by V^, the dressing 
operator associated to the infrared regular fibered Hamiltonian, which has 
therefore the property that 

H''^ = ip2 ^ 1 + V„{x)H,V,{xr + E,{x). 

Using ^e^ai we define the transformed Hamiltonian 

H^^ := %,„H^'''%:, (19) 

which can be expanded in a series of powers of e in J^), with coeffi- 

cients which are also at most logarithmically divergent in a. Thus the gain in 
switching to the representation defined by ^^,(7 is twofold. First we can eas- 
ily separate dressed electrons from real photons and second we can expand 
the Hamiltonian in powers of the small parameter e and thereby separate 
different physical effects according to their order of magnitude. 

The first result we find is that, even though the fibered Hamiltonian has 
no eigenvalues and no spectral gap, there are approximate M -photons dressed 
subspaces which are almost invariant for the dynamics. 

Theorem (Adiabatic invariance of M-photons dressed particles subspaces). 
(see theorem\^ and remark\B\j 

Given any x ^ C*^(M) and any function a{£) such that 

a(£)^/2£"^ ^ 0, e^\og{a{£)-^) 0, e ^ 0+ , (20) 

then 

\\[e-'''''/^,PlMH')\\cijf) < CVMTl\t\e^\og{a{e)-^), (21) 

where 

P'M:=%Ue)i'i-^QM)%Me)- (22) 

The physical mechanism which leads to the almost invariance of the sub- 
spaces is adiabatic decoupling, i. e., the separation of scales for the motion 
of the different parts of the system, which lets the fast degrees of freedom. 
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in our case the photons, instantaneously adjust to the motion of the slow 
degrees of freedom, the electrons. It is a well-known fact however, that the 
decoupling becomes poorer and poorer when the kinetic energies and thus the 
velocities of the heavy particles grow. The quadratic dispersion relation for 
the electrons allows them to become arbitrarily fast, therefore the decoupling 
holds uniformly just on states of bounded kinetic energy. This is the reason 
why we introduce a uniform bound on the total energy of the system through 
the function x, which obviously implies a bound on the kinetic energy of the 
electrons. 

For the following we fix the function a{e) in some way compatible with 
fl20|) . say a{e) = e^. Then we can approximate the dynamics of the particles 
for states in the range of in the following sense. 

Theorem (Effective dynamics of the particles), (see theoremlE) 

Let u be a (mixed) dressed electrons state with finite energy and a fixed 
number of real photons, i. e. u E J^i{Pljx{H'^)J^) and 

uJd:=tr^{%uJ%*) (23) 

the partial trace over the real photons, i. e. the reduced dressed electron den- 
sity matrix. Let the time evolution of uj be the full time evolution 

cj(t) :=e-"^'/"^e'*^'/% (24) 

and define the effective time evolution of Ud by 

uJd{t) :=e-'*^^ff/"cUde"^off/% 

with effective dressed electrons Hamiltonian H^^ given below. 

Let a e C^(M^^ x M.^^) be a "macroscopic" observable on the classical 
phase space of the electrons and Opf^(a) its Weyl quantization acting on 
L2(R3^). Then 

tr^(^(Opf (a)®l^)a;(t)^ = tr^^^jaiv) (^Opf (a) a;d(t)) + 0{e^/^\t\){l - 5mo) 

+ 0{e'\og{e-'m + \t\')), (25) 

where 5 mo = 1? when M = 0, otherwise. 

The effective dressed electrons Hamiltonian is infrared regular and given 
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by 



N 



^,?Jm3~^'^ \k\ 



i,j=i 



2 

with := 1/(1 + y^j) '^"'^ 



er.= ^f . (27) 

Remark 1. The Hamiltonian is equal to the Weyl quantization of the 
Darwin Hamiltonian, which, as we mentioned above, appears in classical 
electrodynamics when one expands the dynamics for the particles including 
term of order e"^ = (f/c)^. At leading order the dressed electrons interact 
through instantaneous pair potentials E{x) given in (USD. At order the 
mass of the electrons is modified (renormalized) and a velocity dependent 
interaction, the so called Darwin interaction, appears. 

Remark 2. The statement fl^^j) remains true if one replaces in ([231) dress- 
ing operator by its leading order approximation V^. 

The dressed electrons subspaces are only approximately invariant 

and transitions between them correspond to emission or absorption of real 
photons. The following theorem describes at leading order the radiated part 
of the wave function for an initial state in the dressed vacuum. 

Theorem (Radiation), (see corollary IE and the subsequent remark) 

For a system starting in the approximate dressed vacuum (M = 0) we 
have, 

^Ut) ■■= (1 - P^)e''^''' P^xm^^ = (28) 
-j^P! K e-ltfV^+^(-)) e-^l^l X: ^p?- 'l ds ^i\'^0^{^{s- X, p)) 

+R{t,e), 
where 

\\Rit,s)\U < CeHogiE-')i\t\ + ItHdl^dll,^' + \\\x\MU + IllPl^dlU) , 
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Xj is the solution to the classical equations of motion 



(29) 

x'jj{0;x,p) = Xj, x''j{0;x,p) = Pj, j = l,...,N, 
K := k/\k\ and 

^Pd:={V,nF,xiH')^)^ (30) 

is the projection of the initial state on the dressed vacuum in the field com- 
ponent. 

Remark 3. Generically (for the precise meaning of this we refer to corollary 
[6]) the norm of the radiated piece is bounded below by 0(^e\og{e~^)), which 
means that the subspace Pq is near optimal, in the sense that the transitions 
are at least of order O (^e \og{e^^)) . 

Remark 4. The radiated wave function \E'rad('^) as given in (PH|) has no limit as 
e because limg^o <^e/|^|^''^ = 0o/\k\^^'^ ^ L^(]R^). However, the radiated 
energy, i. e. the energy of the real photons in (l28l) has a limit. Indeed, if 
we compute the time derivative of the radiated energy in fl28l) we obtain the 
quantum analogue of the Larmor formula from classical electrodynamics. Let 

^rad(t) := (^rad(t), K Hf ^rad(t)) , (31) 



then 

d 3 ^ 

Pr^it) := ^^rad(t) " ^ 5Z ^^^^(^d, Opf (i,^(t; x,p) ■ i^^(t; x,p)) V-d) 

(7/;d,0pf (|d(t;x,p)|2)^d), (32) 



127r 

where d{t; x,p) is the second time derivative of the dipole moment 



N 



d{t;x,p) := ^eiXi{t;x,p) 



i=l 



and the symbol = means that we keep just the leading order in the expansion 
in powers of e. A formal derivation of equation (l32l) is given in remark [8l 

The technique used for proving these theorems is based on space-adiabatic 
perturbation theory, a general scheme designed to expand the dynamics gen- 
erated by a pseudo differential operator with a semiclassical symbol |NeSo| 
|PSTi| ITesI - However, this method exploits the assumption that the fibered 
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Hamiltonian has a spectral gap, of the kind showed in figure [Tj therefore 
is not directly applicable and needs some modifications. In particular, the 
infrared cutoff a plays in our context the role of an effective gap. 

In the more usual framework of time- dependent Hamiltonians, adiabatic 
theorems without gap condition were first proven by Bornemann [Bo] and 
Avron and Elgart |AvEli| (their proof was simplified and the result somewhat 



strengthened in |Tei| ). 

Using similar techniques, a space-adiabatic theorem without gap for the 
Nelson model was proven in |Te2| . In the context of quantum statistical me- 



chanics, an adiabatic theorem without gap for the Liouvillian (the generator 
of the dynamics of the states in a suitable representation of the von Neumann 
algebra associated to the system) is discussed in |AbFri| , while in |AbFr2| a 
general time-adiabatic theorem for resonances is proved. 

All these results are of the first order type, i. e., they describe the leading 
order adiabatic evolution and give upper bounds for the transitions. 

In our case instead, we build adiabatic dynamics including terms of the 
second order, and give an explicit expression for the non-adiabatic transi- 
tions, which allows us to give also a lower bound for them. In the time- 
dependent case, lower bounds for the transitions were calculated by Avron 
and Elgart |AvEl2 for the Friedrichs model, which describes a "small sys- 



tem" , whose Hilbert space is one-dimensional, interacting with a "reservoir" , 
whose Hilbert space is L^(M_,_, (i/i(A;)), with a suitable spectral density /x. 

Finally let us compare our results with those obtained by Davies |Daj . 
He considers the limit c — oo for the Hamiltonian 

:= H^®1 + 1® dV{c\k\) + ^cHi 

which is equal to H^^^ not setting c = 1 as we did before. Davies proves 
then that for alH G M, 

lim e-'^^*(V' ® fir) = (e-^(-^/2+^(-))V) ® • 

c— »oo 

This shows that, while at the classical level the limit £ — > and c ^ oo 
are equivalent, they differ at the quantum one. Indeed, the limit e ^ 
is singular, because no limiting dynamics for £ = exist. Moreover, the 
effective dynamics we get refer to dressed states, which contain a non zero 
number of photons, while the c — oo limit is taken on states which contain 
no photons. 

We summarize the structure of our paper. After some basic facts about 
the model recalled in section [21 we explain how to construct the approximate 
dressing operator ^ in section [SI In section [H we analyze the expansion of 
the transformed Hamiltonian and then apply these results to the study of 
the effective dynamics and the radiation in section [5l 
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2 Description of the model and preliminary 
facts 



In this section we complete the presentation of the Nelson model, discussing 
also the spectrum of the fibered Hamiltonian, and collect some basic facts 
we will use in the following. 



2.1 Fock space and field operator 



(The proofs of the statements we claim can be found in ( |ReSi2| , section 
X.7)). 

We denote by the subspace of the Fock space, defined in ([6]), for 
which ^f(^) = for all but finitely many M. Given / G L^(R'^), one defines 
on ^fin the annihilation operator by 

{a{mY''\h,...,kM):=VMTl [ dkf{k)¥''-^'\k,h,...,kM). (33) 



The adjoint of a{f) is called the creation operator, and its domain contains 
^fin- On this subspace they satisfy the canonical commutation relations 

[aif),aig)] = 0, [a(/)*, a(^?n = . ^ ' 

Since the commutator between a(/) and a(/)* is bounded, it follows that 
a(/) can be extended to a closed operator on the same domain of a(/)*. 
On this domain one defines the Segal field operator 

$(/):=^K/) + «(/)*) (35) 

which is essentially self-adjoint on ^fin- Moreover, ^fin is a set of analytic 
vectors for $(/). From the canonical commutation relations it follows that 

[$(/), $(^7)] = i53(/,(7)i2(M3). (36) 

Given a self-adjoint multiplication operator by the function u on the 
domain D{uj) C L^(R^), we define 

^^,fi„ := C{Qf, a{hy ■ ■ ■ a(/M)*fiF : M G N, G D{u),j = 1, . . . , M}, 

(37) 

where C means "finite linear combinations of". 
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On ^tj.fin we define the second quantization of tu, dr(Lo'), by 

M 

(drHvI/)W(fci, . . . , ku) ■■= J2^{k,)^^''\k,, kM), (38) 

i=i 

which is essentially self-adjoint. In particular, the free field Hamiltonian H{, 
equation acts as 



M 

iH,^Y''\k,, ...,kM) = 5^|fc,|vI/(*^)(A;i, . . . , kM) 



and is self-adjoint on its maximal domain. 

From the previous definitions, given / G D{uj), one gets the commutation 
properties 

mu), aim = a{uf)\ [dr(a;), a(/)] = -a(a;/), 
[dr(.;),i$(/)] = $(i^/). ^ ^ 

2.2 The Nelson model 

Using the Segal field operator (which involves taking a Fourier transform) 
the interaction Hamiltonian, equation ([8]) and ffTOl) . can be written as 

N 

Y,Hi,{x,) = ^{\k\v{x,k)), (40) 
i=i 

where 

„(.,*):=f:e'-.|M.f^e.-.e,|||. (41) 

This form is useful to prove some standard properties of and iffib(a;). 
Lemma 1. 1. is self-adjoint on (see eq. f llll) ). 

2. For every x G M^^, Hfiiy{x) is self-adjoint on D{Hf). 

Proof. The claims are based on the standard estimates (see, e. g., [Be], propo- 
sition 1.3.8) 

||$(/)^||.^<2||//v^||i.(i,3)(V^,i/fV^)^^ + 2||/||i.(K3)||^ (42) 
for G D{Hf), and 

mf)^\\% < 2||//v^||i.(j,3)(vl/, (1 ® i/f)^)..^' + 2||/||i.(M3)||^r^, (43) 
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for \E' G L^(R^^, D(iff)), which imply respectively that, for fixed x, k)) 
is infinitesimally small with respect to H{ and, as an operator on J^, is in- 
finitesimally small with respect to H^, uniformly in e. □ □ 

The fibered Hamiltonian has the form of a quadratic part, which corre- 
sponds to the free field, plus a term linear in the annihilation and creation 
operators. Hamiltonians of this form are usually called in the literature "van 
Hove Hamiltonians" (a review on the subject is given in [De] ) . 

Their simple form allows to analyze in detail their spectral and dynamical 
features. 

In the finite dimensional case, if one considers an harmonic oscillator with 
a linear perturbation, like an external electric field for example, the natural 
way to recover the spectrum of the Hamiltonian is to translate the x variable, 
transforming the initial Hamiltonian into a purely quadratic one. 

In quantum field theory, the analogous strategy would be to find a unitary 
operator which translates the annihilation and creation operators. Such op- 
erator would implement what is called a Bogoliubov transformation. While 
for a finite number of degrees of freedom the Bogoliubov transformation is 
always implementable, this may not be the case if the phase space is infinite 
dimensional. In particular this is not possible if the van Hove Hamiltonian 
is not sufficiently regular in the infrared region. 

After this preliminary remarks, we can state 

Lemma 2. Given the fibered Hamiltonian 

where z{x, k) satisfies 

a{x) := -- [ rfA;^%^ < oo Vx G M^'^ 
2 Jrs \k\ 

one has: 

1. The spectrum of Ksb{,x) is given by [a{x),+oo), and 

(Xac(i^fib(a;)) = (a(x), +oo) . (44) 

2. The infimum of the spectrum, a{x), is an eigenvalue if and only if 

dk — — -- — < oo 

w 
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In this case, moreover, the unitary operator 

V{x) := e'*('^("''=)l'=l"') (45) 
is well-defined on the Fock space and 

irfib(x) = V{x)H{V{xy + a{x). (46) 
If z{x, k)\k\~'^ ^ L^(M^, dk), then the spectrum is absolutely continuous. 



Proof. Point 1 is proposition 3.10 of [Dej and point 2 is proposition 3.13 of 
the same paper. □ □ 

Corollary 1. cr(iffib(a^)) = o"ac(-f^fib(a^)) = [E{x), +oo) , where 

E{x) = -- f dk\k\\v{x,k)\^ . (47) 

A more explicit expression for E[x) is given in f|T3|) . We note also for later 
use that E is a smooth function, hounded with all its derivatives. 

Corollary 2. The infrared regularized fibered Hamiltonian i7fib,o-(a^) (see 
equation ( |T71) ) can be written as 

Hfib,.{x) = V^{x)HfV^{xy + E^{x), 

where 

V^{x) := e'*^'"'^^"''^" . (48) 

Efj{x) is the only eigenvalue, with eigenvector Q„{x) := V^{x)Qf. 
Moreover, for every a G N^'^ , 

d^E, = d^E + 0((tI"I+1)£(.^.) . (49) 



Proof. The only thing left to prove is equation fl49|) . which follows immedi- 
ately from the fact that 

d'^E{x) - d'^E,{x) = [ dk\k\d^\v{x,k)\^ 

2 J\k\<cr 

and 

\k\d^\vix,k)\^ < C|A;|I°I-^ 

□ □ 

As last preliminaries, we show that we can approximate the true dynamics 
with infrared regular ones and that the same holds for the cutoffs in the 
energy. 
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Proposition 1. 



e 



-<=-'"'"'-|L«./^,,<C|(|—. (50) 



e 



norm. 



where M'q = Diiji^^l"^^ with the corresponding graph 

Proof. Both Hamiltonians are self-adjoint on J^o, therefore, given \1/ G Jifo 
we can apply the Duhamel formula to get 

^ Jo 

= - I rfse^(^-*^^'"'/^$(|fc|t;(x,A;)l(o,.)(A;))e-'^^''/^^ 
^ Jo 

^ ||(e-«^V.-e-*^^'^/^)vl/||^< 1 f d.s\\^{\k\v{x,k)\oAk))e-''''''^n.^ 

^ Jo 

where l(o.o-)(^) is the characteristic function of the interval indicated. Using 
equation (HHj) we find that 

||$(|A;|t;(x,A;)l(o,.)(A;))e-^^^/^vl/||^< 

< v^|||fcr/2t;(x,A;)l(o,.)(A;)|U.(M3)||(i/o')^/V'^^^/^v[;||^ 

+ V2\\\k\v{x, k)l(o,a)ik)\\L2(u3)\\^\\jf . 

Moreover, for (3 > 0, 

N 

/c)l(o,a)(/c)||L2(R3) < ^\ej\ ■ |||A;|^~^/V(^)l(o,a)(^)||L2(i;3), 

(51) 

|||A;|^-3/2^(fc)l(o,.)(A:)||i.(^3) = ^^. 

Since ^{\k\v{x, k)) is infinitesimally small with respect both to 1 (8> iff 
and Hq (as it follows again from P3l) ). the graph norm defined by (H^Y^'^ 
and the one defined by {H'^Y^'^ are equivalent, uniformly in e (theorem X.18, 



ReSi2| ). This implies that, for a sufficiently small, 

\\^\k\v{x,k)l^o,.m)e-'''''^'n.^<C^^^^ 
which proves the statement. □ □ 

Lemma 3. Given a function x ^ C'o*^(^^) then 

\\xm-xiH^n\\ciJ^)<Ca'/'. (52) 
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Proof. Using the Hellfer-Sjostrand formula (see, e. g., DiSj| chapter 8), given 
a self-adjoint operator A, we can write 



x{A) = -f dxdydx''{z){A-z)-\ z:=x + iy, (53) 

where ^ C'^(C) is an almost analytic extension of Xi which satisfies the 
properties 

ViVGN : \dx''\<Dfi\'^zf , 

(For the explicit construction of such a see |DiSj| ). 
Applied to our case (1531) yields 



X{H') - xiH^n = - [ dxdy dx^iz) [{H' - zY' - (H^''' - z)''] . 
Since both Hamiltonians are self-adjoint on J^q we have 



{H' - zy^ - {H'^" - zy^ = {H'^" - zy^H'^" - H'){H 



and hence 



Wxm - xiH^nwci.^) < 

< - [ dxd2/|9x"(z)|||(if^''^-z)-l^(^0ll'^(l^l^l(o,.))(^'-^)"'IU(^)- 
In addition we have that 



C 



||(/f^''^-2;)-i|U(.^.)<-— . (54) 

This follows because Jifo is dense in the domain of i/^=°''"=° = L'^{R^^, D{Hi)), 
and for every \I' G 

H^,^^ i7^=o,'^=o^ as (£, a) (0, 0) . 
According to theorem VIII. 25 |ReSii| , this implies that 

therefore I'^zlW {H'^''^ — z)~^'$\\ is bounded for every and the uniform bound- 
edness principle gives ([5] 
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For the second norm we find that for z G supp 

||$(|A;|^;l(o,.))(//^-^)-i£(^) < 

< ||<l>(|A;|t;l(o,.))|U(,#6,..^o-|l(^'--2'"' 



< -^\\^{\k\v\ 



I 1 



io,^))\\c{.jro,je), 



because if^ and define the same graph norm uniformly in e, as it follows 
from ( |43|) . From the same equation and from ( ISTl) it follows also that 

||$(|fc|t;l(o,<x))IU(,«,^) < |||A;r^^^^l(o,<7)||L2(R3) + y2|||A;|fl(o,a)||L2(R3) < Da^^^, 

which proves the statement. □ □ 



3 The approximate dressing operator 

In this section we construct the approximate dressing operator ^e^o- for the 
infrared regularized Hamiltonian H'^''^. 

We recall at the beginning the formal procedure to build the so called 
superadiahatic projectors, introduced first by Berry [Berj for time- dependent 
Hamiltonians and generalized by Nenciu and Sordoni |NeSoj to the space- 
adiabatic setting. 

The basic idea is to construct a sequence of projectors -P/, j € N, whose 
commutator with the Hamiltonian is of order 0{e^^^). Then the range of P? 
is almost invariant for the corresponding time evolution in the sense that 

[e-'*^'/^P/] = 0(£^>|). 

As we have briefly mentioned in the introduction, this procedure, which 
allows to build adiabatic dynamics beyond the leading order, has been imple- 
mented till now only for fibered Hamiltonians whose spectrum is made up of 
eigenvalues separated by a gap. Exploiting the fact that the infrared cutoff 
a gives rise effectively to a gap, we can however carry out the calculations 
also in our case. 

We proceed first in a formal way, using the superadiabatic projectors just 
as a formal tool to deduce the expression of the unitary intertwiner We 
will then show that the expression we find for ^ actually defines a unitary 
operator on ^ and we will prove some of its useful properties, which allow us 

""^To simplify the notation, we drop from now on the dependence of the unitary on a 
and e, denoting it simply by instead that by '^e.ai as we did in the introduction. 
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to fully characterize the domain of the effective Hamiltonian, and to expand 
it in a series of powers of e which is convergent in C{^q, M'\ 

The projectors Pm on the almost invariant subspaces are then recovered 
from ^ via Pm = '^*Qm'^- 

It is important to stress that, to characterize the transitions between the 
different almost invariant subspaces, we need a unitary which maps all of 
them simultaneously to the corresponding reference subspaces. For this rea- 
son, we cannot use the procedure based on Nagy formula proposed in |MaSoj 
|NeSoj and we employ instead a simpler one, based on the construction of 
the projectors. 

In what follows, a is considered a parameter independent of e. Only at 
the end, we will choose it as a suitable function of e. 



3.1 Formal definition 

Using the dressing operator Vfj{x), whose expression is given in (HHI) . we can 
define dressed M-particles projections: 

7r*^(x) := V„ix)QMV^{xr. (55) 

Note that, for M = 0, we have Il^\x) = |fio-(a;))(^^o-(x)|, the projection 
onto the ground state of Hfi^^^^^x), but, for M 7^ 0, (x) is not a spectral 
projection of the fibered Hamiltonian. 

The standard construction we briefly mentioned above, which is used in 
the case with gap, is based on an iterative procedure. Starting from a zero 
order projection n*^") = ttq, which corresponds to a spectral subspace of the 
fibered Hamiltonian, it allows to build an approximate orthogonal projection 
n*^"^ := X]j'=i^"'^j which satisfies 

[^.^n(")] = ^^^^ 

As we already stressed, we avoid on purpose to be more precise on the 
sense in which equation ( l56i) holds and on which norm one should use on the 
right-hand side, because we will use it just as a formal tool to deduce the 
expression of 

To clarify the method we employ, we recall briefly the derivation of the 
first order almost projection II*-^-* in the case with gap, i. e., we assume for 
the moment that Hsh^^{x) has a nondegenerate eigenvalue Ea^x) which is 
isolated from the rest of the spectrum. Our discussion is taken from IPST2 . 
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Therefore, we start from a projection tTqIx) which projects onto the 
eigenspace associated to E„{x), and we proceed formally, without worrying 
about specifying any regularity assumption. It holds then that 

K = ^0. (^o)^ = (57) 

[tto, = [^o> I P'] = i^lViTo) ■ p + y Atto, (58) 

where the right-hand side of the commutator is of order 0{e) when applied 
to functions of bounded kinetic energy. Applying the inductive scheme, we 
determine now the coefficient tt^ in II*^^-* = + ^^i by requiring that 

(n«)* = n«, (n«)2-n« = o(£2), 

[n«,if^] = o(e2). ^^^^ 

The first condition gives 
so that we must have 

TT^ = TTpVT^ + Tr^TTg + 0(e). (60) 

Concerning the commutator, we have 

[U^'\H''''] = ie{Vno) ■ P + e[n„ i/fib,.] + 0{e') = 0{e'), 
so has to satisfy 

[tti, i7fib,a] = -i(V7ro) • p + 0(e). (61) 

It can be shown ( [Teal , lemma 3.8) that the unique solution, up to terms of 
order 0{e), of (1611) is given by 

TTi = -i7ro(V7ro)(i/fib,. - Ey'n^ ■ p + adj. , (62) 

where ttq := 1 — tiq, and + adj. means that the adjoint of everything which 
lies to the left is added. This ensures also that tt^^ is formally self-adjoint. 
If there is a gap, the reduced resolvent is certainly bounded in norm: 

(63) 



inf dist ( E^{x), a{H{n,^aix))\{E^{x)} 
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but in general the construction breaks down in the case without gap, because 
the reduced resolvent can be unbounded. 



A possible way out of this difficulty was proposed in [Teij and |Te2| to 
prove the adiabatic decoupling for the zero order subspaces. It consists in 
introducing a second parameter (5 > 0, by which the reduced resolvent is 
shifted into the complex plane: instead of (-fffib — E)~^7r^, one considers 

In our case, the infrared cutoff a plays an analogous role. Actually, at 
the end of our formal calculations, it will result that the reduced resolvent 
of -fffib.o- is bounded. 

We carry out another generalization of the standard construction, extend- 
ing it to the situation when M 7^ 0. To understand the way how we proceed 
it is useful to analyze more explicitly the case M = 0. According to the 
above discussion we define 

n«(e,a) := IlP+en',{a), 

7r?(a) := -i(Vn(°))i?4(E.) • p + adj. , 



where R^^{Ea-) := (-f^fibo- — -Eo-) Hq . We have omitted the IIq because 



of the well-known fact that Vfll'^'' is off-diagonal with respect to the block 



decomposition induced by nQ^"*. 

Using the fact that VnJ,"^ = iVa{x)[^{iVa:V^) , Qo]Vcr{x)*, and introducing 
the notation -Rj^(O) := H^^Qq, we get then 

= -V^ix) [dr(|A;|-i), [$(iV,t;.), Qo]] ■ pVAx)* + 0(e), 

where we have used that the commutator between p and a function of x is 
of order 0{e) (we stress again that we don't worry here about smoothness 
assumptions, the calculations should be considered as formal) and 

RiiO)Q<i = dT{\k\-^)Q<i, 

together with 

Qo*(iV^^^a) = QoHiVxVa)Q<u ^{iV^v^)Qo = (5<i$(iVxWa)Qo- 
Calculating the commutator using equation fl39|) . we get 

n',=iV4Qo,^^]-pV:, 

where we define 

<bJ'r^ ■= (bf^ 

\k\ 



$.(a:):=$(^^^#^'). (64) 
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Since V xVa{x, k)\k\~^ G L"^ (M^^ , dk) , ttJ' is well-defined, and one can sym- 
metrize it to get a symmetric operator (the fact that 7r° is not bounded due 
to the presence of p will be dealt with below). 

It is now fairly clear how to proceed. For any M, we define 

i (65) 

Formally II^'* satisfies (!59|) for any M. In fact, $0- is self-adjoint, therefore 
i[QM,$o-] is self-adjoint, and (H^lj)* = Moreover, since vrf^ is off- 

diagonal with respect to tTq^, 

^0 ^1 ^0 = (1 - ^0 (1-^0 ) = 0' 

equation (I6OII holds exactly, without 0{e) corrections. 
Concerning the commutator equation ( 16T1) . we get 



[yrf , i/fib,,,] = [[Qm, i$.] ■ P, ^^f + ^(x)] + 0{e) = 

= V^[[QM,i^a],Hi]] ■pV: + Oie) = V^[[Hi,i^^],QM] ■pV: + OiE). 

Applying again equation flHIH) we get, on ,^\k\fin, 

[iff,<l>.] = <l>(iV.t;.), (66) 

so that 

[vrf , i/fib,.] = -iVTT*^ ■ p + Oie). (67) 

The next step is to find a unitary operator which intertwines the almost 
projections II^'* with the reference projections given by the Qm up to terms 
of order ^(e^). 

Employing a procedure analogous to the one we used for the projections, 
we assume that we can write an expansion 

n 

U(^) ■=J2e'Uj, (68) 
i=i 

starting from a known Uq, and imposing that 
to deduce the coefficients [/„ iteratively. 
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An obvious choice for the zero order unitary is given by 

f/o(x) := (70) 

which satisfies Uo{x)Il^^ {x)Uo{x)* = Qm- 

Without loss of generahty, we assume that Ui = TUq, for some operator 
T on J^. The requirements ( !69l) then lead to 

{Uo + eUi){U* + eU*) = 1 + e{UoU* + U^U*) + 0{e^) = 

= l + e{T* + T) + 0{e^) = 0{s^) =^ T* + T = 0{e), 

and 

(Uo + £t/i)(n£) + e7rf){U* + eU*) = (1 + eT){QM + bIQm, i<f .] ■ - eT) 
+ Ois^) = Qm + e{[QM, ■ P - [Qm, T]) + 0{e'') . 

Therefore, by choosing 

T = i$.-p, 

we satisfy both requirements for every M. The first order almost unitary is 
then given by 

= (i + i£$<,-p)r;. (71) 

3.2 Rigorous definition and properties 

To give a meaning to the till now formal expression for U^^^ we have to deal 
with two problems. 

The first is due to the fact that the operator $o- is unbounded. We intro- 
duce therefore in its definition a cutoff in the number of particles, replacing 
it by 

:= Q<j$<.Q<j. (72) 

The operator 

f/f := (l + ie$^p)K* 

is again formally almost unitary up to order O(e^), and intertwines the su- 
peradiabatic almost projectors for M < J, 

t/f n£)f/f* = QM + 0(e2), M<J. (73) 

This means that we can use f/j^"* to study the transitions among superadia- 
batic subspaces up to an arbitrary, but fixed J. 

The second problem, already mentioned in the introduction, is related to 
the presence of the unbounded momentum operator p. 
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To make the whole expression bounded we introduce a cutoff in the total 
energy. More precisely, given a function x ^ C*^(M), we define 

ufi := v:[i + iev^^i ■ pv: - i£(i - xiH'nWa^i ■ pv:ii - xiH^n)] = 
= v: [1 + iexiH^nvM ■ pv: + hi - xiH^nw^K ■ pvMH^nm) 

Note that U^^lxiH''") = U^'^^xiH''"), for every x e C^{M.) such that XX = 
X- In the context of Born-Oppenheimer approximation, Sordoni |So] applied 
a similar strategy in order to define bounded superadiabatic projections. 

Our aim in this section is first to prove that Uj^^ e C{J^) nC{Mo). Once 
we have shown this, we can define a true unitary through the formula 

We will then prove that both and = '^* belong to n /:(^), 

i. e., that ^ is a bijection on We will in addition show that it can be 
expanded in a convergent power series both in and in L^-Wq). 

Lemma 4. K- and belong to C{J^o). Moreover 

\K\\ciM)<C, (75) 

with a constant C < oo independent of a. An analogous estimate holds for 
V* 

Proof. {We give the proof for V^j. The one for is the same up to some 
changes in the signs). 

We can calculate the norm on a dense subset. It is known ( |ReSii| , 
theorem VIII. 33) that Hq is essentially self-adjoint on a set of the form 

Vp ® Vi := C{ij : eVp,^ e Vi}, 

where C means "finite linear combinations of" , Vp is a core for p^ and T>{ is 
a core for H{. We can choose then Vp = C^{R^^) and ( IReSial , section X.7) 

A = {V? G ^fin : v'^'') e ^%,C^{R') n ^^(M^^)}, (76) 

where ®%^C^iR^) := C{ipi (g) . . . (g) cfN : e C^(\^ ^),J = 1, . . . , N}. 

The vectors in are analytic vectors for V„{x) ( |ReSi2| , theorem X.41), 
so, if ip E Vp, G V{, we have 
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Moreover 



- 2iej^w^y-^^{iV^v^) ■ pil){x)^ + ^{}v„yp^i){x)^ , 

where we have used (IHUj) and the fact that '^{v„{x, ■), Va;fo-(a:, ■))l2(r3) = 0. 
From the previous equations it follows that 

is convergent, so that Va^ip ® ip) C DijP' ® 1) and 

(p2 ® i)v;(^ ® v^) = -i£V,$(iAi;,)(V' ® (^) + e^V^^iiV^v^y (77) 
•$(iV^w,)(V' ® (^) + 2eK$(iVx.t;a) ■ P{iIj ® ^) + V^{p^ ® 1)(^ ® (^) . 

This implies 

II (p2 ® ® v^)|| < e^lMiAv^y^ ® <^)|| + 

+ £l<l'(iV,t;^) ■$(iV,t;,)(V^®^)|| +2£||<l>(iV,t;,) ■p(^®(/^)|| (78) 
+ \\{p^0l)iij0if)\\. 

Using equation (H3l) we can bound the first term with some constant inde- 
pendent of 0" and e times ||^/' ® v^|| js^fj. 

This happens because the constants involved in fHHj) contain terms of the 
form ||//iy[fc[||2,2(iR3), but when one differentiates with respect to x, one 
gets an additional |fc|, therefore all the terms of the form || V^fo-/ vT^II 
so on are uniformly bounded in a. The same reasoning applies also to all the 
estimates which follow. 

For the second term in fl77|l . using again fHHj) . and the notation = ip^ip 
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and / = idjV(j, we get 



+m\\hm\mmr^ 

< v^ii//^iii.(M3)($(/)^, Kl ■ IfT - «(l ■ 1/)]^).^ 

VnilL2(R3) ■ ||$(/)\I/||.if' ■ II [a{\ ■ I/)* - a(| ■ |/)]\I>||.^ 
+ 2||//v1^|li.(M3)-||$(/)^^||.^H|i^f^lU + 2||/||V)||$(/)^||^ 

< ^ll//v1^lli2(M3)||$(/)^lL^ + ^11 Kl ■ \fr - «(l ■ 

+ &'ll//v1^lli.(M3)ll$(/)^^ll.^^ + ^ll//v1^lli.(M3)||i/f^l^ 

+ 2||/||i.(M3)||$(/)*||!^, Va,6>0. 



mfrn%il-b^f/^A^m< 
< ^\\f/^Al\\lH^^)mm\\% + J^w Kl ■ ifT - «(i ■ i/)]*ii' 

+ ^ll//v^lliwll^f^ll'^^ + 2|l/lliwll'^(/)^ll^ • 



Each term on the right-hand side can be bounded by a constant independent 

of a times ||\l/||^o- 

For the third term in fl78p we get 



||$(m-*||^ < 2||//v1^||i.(^3)(p,^,i?fP,^) + 2||/||i.(j,3)||p,^||^ < 

< 2||//v^||i.(^3)(vP,i/fP^'v^) + 2||/||i.(i,3)||p,vl/f^^.< 

< ll//v1^lliwll(^l + ^f)^ll^5^ + 2|l/llW^)llp.^lL^. 



and we can again bound the right-hand side by a constant times ||\l/||.^. 

We have now to prove similar estimates for || (1 (8> Hf)Va{ip ® v^) || , and we 
are done. For this we need that, on Vp(^Vf, 
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Hence 



[Hi,^w^)] = imk\v,) 



and that 



We have then 

i-i 

(1 ® Hi)<l>{iv^{x, ■)y{^ ^^) = iJ2 <l>{iv^y<l>{\k\v^)<l>{Ky-'-\^ ® ^) + 

1=0 

= ij^w^y-^mk\v^){ilj ^ip) + E^{x)j{j - l)^{iv^y-^{ilj ® if) 

+ ^{w,y{l(»Hi){ilj®ip). (79) 

This means, as for p"^ 1, that Va{ip ® (p) C D{1 iff), and that 

{l(^Hi)V^{ilj^ip) = ~V^^{\k\v^){ilj®ip)-E^{x)V^{ip0ip) + V^{l®Hi){ilj^ip). 

(80) 

Using equation (143|) . we can again bound each term by a constant indepen- 
dent of cr times Wip ® v^H.^^. □ □ 

Lemma 5. i. For eac/i x G and / = !,..., 3A^, the components 
o/ the operator $^ satisfy ^i^i{x) G £(^) and $^_,(x)* = $^/(x). 

Moreover, 

<I>J,, : M^^ ^ £(^), x^$J,(x), GC^(M^^/:(^)), (81) 
and, for a small enough, 

WKihiM^) < Cv/jTTv/log(a-i). (82) 
Given a G wt/i |a| > 0, /ioWs instead 

d^^i^i = d:^i, + 0(al"lv/jTT)^(^), (83) 

where 

:= (9:$J,,)|.=o (84) 
zs a well-defined bounded operator on Jf. 

2. The statements of point 1 (except for the self-adjointness o/$^;(x)j 
remain true if ^ is replaced by D{Hf). 

Proof. The proof follows applying the standard inequality 

\\Q<Mf{x,kyQ<,\\ci.,r) <2'/'VJTl- sup ||/(x,A;)|U2(K3.,fc), 

to the case / = d"diVa □ □ 
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Corollary 3. The fibered operators 9"$;^; belong to C(J^) fl C{J^) Va G 
Lemma 6. Uj^^ belongs to C{J^) fl C{J^fo) and 



f/SlU(,^) < C(l + £Vlog(a-i)), (85) 



where Ji^ = Jf or J^q. The same holds for Uj^^* . 

Proof. We start considering Uj^^. We have already shown that and V* 
belong to C{Jif) fl C{J^) with uniformly bounded norms, so we have to 
examine 

exiH^nv.K ■ pv:, e{i - xiH^nwK ■ PVMH^n ■ 

Given \1' G Vp V{ (for the definitions, see the proof of lemma H]) we have 

Since V^; ■ (K-$;^) and x(/7^''^)p are bounded operators (see lemma H] and [5]) 
the left-hand side belongs to C{J^) fl £(=^). Moreover, from the previous 
lemma it follows that 



while all the other terms have uniformly bounded norms. 
For the second one, 

v^^i ■ pv:x{H'n'^ = v^^i ■ [p, vMH^n^^ + ■ v:pxiH'n'^. 

px{H^) G C{J^) n C{Jifo), and so do all the other terms, as it follows from 
lemma H] and lemma [5l 

Concerning f/j , from equation ( TMI) it follows that 

U^^l* = [l-iEx{H'nV.p-Kv:-ie{l-x{H'n)V.p-Kv:xiH'n]Va, (86) 

so one can apply the same reasoning as above. □ □ 

Theorem 1. Assume that a = (j{e) and that conditions (J20l) hold, then the 
operator 

^:=<[<*<]-/^ (87) 

is well-defined and unitary, for e small enough. 

Both and belong to C{J^) fl C{J^o), with the property that 

||^IU(,^6)> ll^*ll£(.i^o) <C, (88) 

where C is independent of e and a. 

Moreover we can expand them in a power series which converges both in 
C{Jf) and %n £(^). 
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Proof. Combining equations (17^ and flHBl) . and using the fact that p ■ $^ = 
■ p — ieV ■ we have that 



where 

-Be := Bq + eBi, 

Bo := -x(i^^''^)KV • <^iV: - (1 - x(^'''^))KV ■ 'i>iV:x{H'n 

+ [xiH^nv.iK ■ p)v: + (1 - x{H^n)v.{K ■ pwMH^nnm 

B, := -i[x{H'nV.V-<!>iV: + {l-x{H'n)V.V-<!>iV:x{H'n]- 

■IxiH^nv.iK ■ p)v: + (1 - . 

^From lemma E] it follows that i?e is a bounded operator both on Jif and 
self-adjoint on J^. From equation (!85|) we have also that 

\\Be\\ci.x) = 0{\og{a-')), = ^,^0, 

so, under conditions ( 120|) . we can assume that ||e^Se||£(^) < 1. Uj^^*Uj^^ is 
therefore strictly positive, and the square root in (187|) is well-defined. 

Moreover, we can express it through a convergent power series, both in 
C{J^f) and in 

= il + e^B,)-y- = f:i-iy^^e^^BlE. (90) 
We can also explicitly calculate that 

where we have used in the first equation the fact that both Uj^^ and Uj^^ * are 

invertible, since they differ by a term of order 0{e>y\og{a^^)) from a unitary 
operator. 

Putting together (ISTIl and we get in the end the expansion for . 

□ □ 

i(2j-l)!! = 1.3-5---(2j-l); (2j)!! = 2 • 4 • 6 • • • 2j = 2^j! 
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4 The dressed Hamiltonian 



Using the results of the last section, we can define the dressed Hamiltonian 
just as the unitary transform of H^'" , 

Hl^^^ := ^H^^^^\ (91) 

It follows from theorem [1] that is a bijection on so H^^^^ is self-adjoint 
onJ^. Moreover, since can be expanded in a power series in we 
can expand i^^res -^(=^5 -^)- However, putting directly the expansion for 
in (|9T1) we get e— dependent coefficients, because of the e— dependence in 
H^'" . To get the correct expansion we need then to rearrange some terms, 
but the remainder we get will in the end be bounded in J^). 

Theorem 2. The expansion up to the second order of the dressed Hamilto- 
nian is given by 

Hl^^ = /io + ^/^i,x + ^'/^2,x + (^M£(.if6.ir) + 0(£'(log(cr-i))'/')^(.^^,^0' (92) 
where hi^^ G C{J^q, J^) , 

ho = ^p^ + H{ + Eix), (93) 

^i,x given in equation ( l95l) and /i2,x is given in equation ( 1971) . 

Remark 5. The coefficients in the expansion of -ffdres depend explicitly on the 
cutoff function x (and also on the cutoff J in the number of particles, even 
though we have not stressed it in the notation). 

However, as we have already mentioned, we expect the adiabatic decou- 
pling to be meaningful only on states in the range of xiH'^'") (or x(if^), 
according to lemma (El) we can interchange the two), which becomes in the 
representation space x(-f^dres)- 

With this in mind, we will prove later that the effective dynamics on the 
range of Qa/ x(-f^dres) < J+1) is generated by an Hamiltonian independent 
of X and J. 

Proof. To simplify the following reasoning we write 
Ufl = V:{l + eT), 

T := ix{H'nv.K-pv:+i{i-x{H'n)v.^i-pv:x{H'n. 
s := ix{H'nvaV-<^iv:+i{i-x{H'n)VaV-<i>iv:x{H'n- 
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We omit the dependence on a and £ in T and 5* to streamline the presentation. 
We can then write 

B, = iS-T'^ + ieST. 
Expanding the square root according to formula (!90i) we get 



Hl^^^ = ^H'^"^* = V;*(l + eT) 



l-i-(iS-T^) 



2 2 



We examine now separately the terms coming from different powers of e. 

1, 



(0) = 



V,{xyp'V,{x) + H, + E,{x) 
^p^ + Hi + E,{x) + \v,{xr[p\ K(x)] 



(94) 



2 2 

+ y$(iV^t;^) • $(iV^w^) - iy$(iAt;, 



ggj 1 



e 



e 



p' + Hf + E{x) + -$(iV^.w^) -p+'P- $(iV^w,) 



This implies that 



Summing the term of order 0{e) coming from (0) to the ones coming 
from the expansion of ^ we get 

(1) = ^^(iV,^;.) -p+^-p- H^.v,) + eV:[T, H'''^]V, . 
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Here 

[T, H''"] = 

= '-x{H^nK^i-pv:,p'] + ^(1 - x{H^nw.<i>i-pv:,p'MH^n + 
+ i(i - xiH^nmK ■P,Hf+ E^{x)]v:x{H'n ■ 

The commutator appearing in the last two hnes gives 

[$J ■p,Ht + E,{x)] = [$J , ■ p + $J ■ [p, E,{x)] = 

= iQ<j$(iV,t;,)Q<j ■ p - ie^i ■ VE„{x) , 

where we have used equation (|66l) to calculate the first commutator. 
Putting together all the terms we have therefore 

(1) = |$(iV.t;.)-p+|p-$(iV.t;.)-£x(K*^''"K.)[$'(iV.t;.)-p 

~e^i ■ VE,{x)] - e{l - x{V:H''''V,))[^' (iV ,v,) ■ p 

-e^i ■ VE,{x)]x{V:H'^''V,) + 0{e^^\og{a-^))ciM,^) , 

where the terms containing the commutator with p^ give rise to higher order 
contributions and we abbreviated := Q<j$(/)Q<j- 

Analyzing the terms of order 0{e'^) we will see that they yield, as it 
happens for the zero order one, terms of the form p ■ $, which make the 
previous expression a symmetric operator. We have therefore in the end that 

/zi,^ = ^$(iV.t;,)-p-^x(K;^''"K)$''(iV.t;.)-p (95) 

- i(i - xiy:H'^''v,))^'{iv^v,) ■ px{y:H'''^v,) + ■ $" , 

where the symbol "p • $" means that, associated to each term of the form 
$ ■ p, there is another one of the form p ■ $ which makes the sum a symmetric 
operator. 

To calculate /i2 we follow the same route, 
(2) = y$(iV.t;,) • $(iV.t;.) + £\(V;*^''"K)<f ' ■ VE„{x) 
+ e\l - xiV:H''''V,))^i ■ WE,{x)x{V:H''''V,) + 

■[^i ■ p, v:p'v,] + |(i - xiv:H^'''v^))[^i ■ p, v:pWmv:h^^''v^) 
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We examine separately the different terms. 

2 2 

^2 3iV 

= jYl H^^'^)'^ + + a{idjv^)*a{idjv^)] + || Vtv|li2(K3,rffc)^c3iv • 

We remark that the last term, being x-independent, is a c-number. 

■ P, v;y K] = [^i ■ p, p'] + [<i>^ ■ p, v: [p\ K]] = 

= 2ied^<^i^ipipj + e[<l>i ■ p, $(iV,t;^) -p + p- ^iV^v^)] + 0{e^)c(M,.^) 
= "^^ed^KiPiPJ + ^iKi^ Hid^v)]pipj + epipj[^i^i, Hidjv^)] + 0{e^)ciM,.^) 
= 2i£9^.$i;PiP, + e[^i^i, ^id^v^)]PiP, + ePiPj[<!>i^i, ^id^v^)] + C(£')£(^o,^) 
+ 0{ea)ci.M>,^), (96) 

where we have used lemma [5l 



1 1 
= -^Qja{id^v^)<^iQ<j - -^a{id^v^yQj<^iQ<j 

+ -^Q<j<l>iQja{id^v,) - ■^Q<j<i>ia{id^v„yQj 

+ iQ<j^d^V^ix, ■), ' " ' )z.2(M3,,fc) 



=: ^j + ig<jg?(g,t;.(x,-), ^'''|jp"0 

/^j + iy<jJK( |_|i/2 ' |.|l/2 /L2(M3,dfc) + ' 

where we have bounded the scalar product using the same procedure applied 
in lemma [5l 

The remainder term vanishes when applied to states in the range of 
Qm with M < J — 1. The scalar product can be written in a clearer way using 
the explicit expression of the function v {h,ji = 1, . . . N, l2,j2 = 1, ... 3): 

\ |.|l/2 ' |.|l/2 /L2(R^dfc) 

= 3? / dk gi-^(k) Qj^{k) j— -- Hi^Kj^, k, := —— . 
Jr3 \k\^ \k\ 
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When li = ji, since the charge densities are spherically symmetric, we get 

1 /■ , , ^Zi(a;)^zi(t/) . 1 . 

J-^TT Jk3xm3 |a; — 1/| 6 

Putting all the terms together, we have in the end 

AT 

aiT^ii) 

where TZj is a remainder term that vanishes when applied to states in the 
range of Qm, with M < J — 1, and Pi_^ denotes the three-dimensional mo- 
mentum operator associated to each particle. 

• 2 2 

2 2 
2 2 

■[iff, V • ^i]x{V:H^^''V,) + 0{e^)ci^,,^) 

• 2 -2 

= ^x(K^^^'"^.)^'(iA^;.) + ^(1 - x{V:H^^''V,))^'{i/^vMy:H'''V,) 
■ 

This gives exactly the terms needed to make hi a symmetric operator. 
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Using the calculations for the first order part hi, we get that 

[//^■^T] = xiH'nva^'i^xVa) ■pv:+ 

+ (1 - x{H'n)Va^'{iV,v,) ■ pV:x{H'n + 0{e^/log{a-')))ciM,^), 
we can keep therefore just the first two terms, omitting the remainder. 

2 2 

+ (1 - x{H'n)ya^'i}^xV,) ■pV:x{H^n.T]V, + 0{e^^\og{a-^)) 

= '^v:[x{H^'nv.^'i}V,v,)-pv:, x{H^'nv.K-pv:]Va 
+ '^v:[{i - xiH^nWa^'c^v^v,) ■pv;x{H'n. xiH^nvaK ■ py:]ya 
+ [x{H'nya^'i:^^:cVa) ■ pv;, (1 - xiH^nKK ■ pv:x{H'n] Va 

+ ^y; [(1 - x{H'n)Va<^\^.v,) ■ pVMH^n, 

(1 ~ xiH'-nwM ■ pv:x{H'n] ^ + v^m^)£(^o,^) • 
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Summing up, we get in the end that 

^ 3N 

h2,x = 4 [aiidjV^f + aiid^v^)*"^ + a{idjV„ya{idjV„)\ (97) 
i=i 

+ l|Vt;.||i.(M3,,fc)«c3" + x{V:H^'''V,)i^^i ■ VE{x) 

1 ^ 

h=i 

'iji=i 

+ («:•pJ(/^•p,.Je^'^■(-^l"-'l)]+7^J| 

+ (1 - xlK^if^-'^K))! • • ■ 

+ ^r;[x(^^''^)v;$^(iv.t;.) -pv;*, x(i^^''^)K<^>^pv;*]v; 

+ ^v: [{1 - x{H^-n)Va^'{:iv^v,) ■ pYMH^n, xiH^nv.K ■ pv:] k 

+ ^v:[xiH'nva<^\iv,v^)-pv:, (i - x{H'n)VaK-pv:xiH'n]Va 

+ V; [(1 - xiH^nWa^'iiv^^v^) ■ pvMH^n, 
(1 - xiH^nWaK ■ pvMH^n] K , 

where we have symmetrized dj^Q^ipipj, which is possible up to terms of order 
0{e). The expression is fairly lengthy, but we will show below that many 
terms vanish when applied to a state in the range of Qmx{H^^^^) (M < 
J-1). □ □ 



5 The effective dynamics 

We start with a number of lemmas we need to analyze the effective time 
evolution. 

Lemma 7. Assume that a satisfies conditions (1201) . then 
1. Given a function x £ C'o°(IR); we have 

x{y:H''''V,) - x{ho) = £7^^, (98) 
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where 7^^ G C{J^,J^o), \\'R-'1,\\c{.j^,M) — ^(1) 

n^QM = (Qm+1 + QM-i)niQM + 0{e^)c(.je,M)- (99) 
2. Moreover, we have that 

X(^drcs) - x(/^o) = (100) 

QMx(i^d'reJ = QMXihomHLs) + 0(5^ V^MO )£(«) , (101) 

where x is any C^(R) function such that XX = X ^^'^ XX = X; M < J — 1. 
Proof. Applying the Hellfer-Sjostrand formula, equation we get 

x(v;*i7^'-v;) - xiho) = - [ dxdy Briz) [{v:h''''v^ - z)-' - (ho - z)-'] . 

Since both Hamiltonians are self-adjoint on Jifo, we get 
{V:H''''V^-z)-'-{ho-z)-' = 

= {ho - z)-\ho - v:H'^"v^){y:H''''v„ - zY^ 



-(/lo - zY 



+ -$(iV,t;,) ■ $(iV,t;<,) 



£(jr,jro) , 



where we have used equation to calculate the difference of the two 
Hamiltonians and estimated the integrand proceeding in the same way as 
in lemma El Moreover, iterating the formula we get 



{y*H^,'^V,-zr' -{ho-z)~^ = 

= (ho - zr\ho - ic^''"K)(v;i/^''^K - z)-' 



(102) 



-(ho-z)-' 
+ 0{e^\'^z\-^ 



{ho - zY' 



so 



',M 



[{V:H'^''V„-z)-^-{ho-z)-^]Qi 

= (Qm+1 + QAf-i) [{V:H''''V„ - z)-^ - {ho - z)-'] Qm 



+ 0{e'\Qz\-')ci 
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Concerning point 2, equation fllOOl) follows immediately from the fact that 

^dres - - (ho - = {ho - z)~\ho -H^rJiHlres ' ^Y' = 

= e{ho - z)-'h,^{Hl^^^ - z)-' + 0{e',/I^i(^)ci.M,.^), 
while, for equation f llOip . we have by definition that 

= QMX{ho)x{Hl,J + Qm[x{HI,J - x{hoM{H'^re.) ■ 
Proceeding as above we find 

If- r 1 n~ (103) 

= - / dxdy dr{z)QM[m,,, - z)-' - {ho - z)-']x{m,j , 

so, if we show that 



QM/ii,xX(^drcJ = 0{e^\og{a-^))ci.^,,,^) (104) 

we are done. 

From equation (omitting the "p ■ $" part, which can be treated in 
the same way) we get 

QM/^l,X^(^drcs) = \QMnV.V„)-PmirJ " ^Qa/ X(K*^'''^K) " 

■ ^\iv.,v„) ■ M(^drcs) - \Qm{-^ - x{v:h''-v,))^'{;iV^v,) ■ Px{v:h''''v,)- 

but it follows from point 1, that we can replace x(y*H'^V) with x(^o) up to 
terms of order 0{e), and from equation fll03p that we can replace x(-f^dres) 
with x{ho) up to terms of order 0{e), therefore we have 

QA//il,x^(^dres) = .V„) ' Pmo) " ^QmX(K*^''"K) (l V.^;.) ■ 

mho) - IQm{1 - x{V:H''''V^))<^'{W,v^) ■ p^{ho) + 0{e) 
= ^QaMiV^v^) ■ P^{ho) - ]^QM^'iiV.v^) ■ P^{ho) + 0{e) = 0{e) , 
so point 2 is proved. □ □ 
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The following lemma was proved in ( [FGSj . Appendix B) and we will use 
it to characterize the range of X(-oo,c)(^o)) where X(-oo,c) is the characteristic 
function of the indicated interval. 

Lemma 8. Let H be a Hamiltonian of the form 

H ■=l®dT{\k\) + H®l, (105) 

acting on the Hilhert space M' = 3^ ® ^ , where ^ is the bosonic Fock space 
over L^(M'^) and is a generic Hilbert space. 

Then the set of all linear combinations of vectors of the form 

N 

<^®a{giy ■■■aig^yriF, A + ^Mj<c, (c> 0), (106) 

i=i 

where ip = X{-oo,\){H)f for some A < c, N G N and 

Mj := sup{|A;| : k G supp gj} 
is dense in X{-oo,c){H)J^ . 

Lemma 9. Let x ^ C^i^) andho := 1® H{ +p'^ /2®1. Then, 3^ G C^(K) 
such that xi = X 

nh,)x{h,) = 0{e^)ci,^), (107) 

where ^'^ := 1 — S,- 

Moreover, denoting by 

:= sup{|A;| : k G supp 

and defining 

:= 2c^ + Eoo, (108) 

where Eoo '■= sup^g]R3n|-E'(x)|, we can choose sup{|fc| : k G supp ^} arbitrarily 
close to cg. 

The statement remains true also if we invert the roles of Hq and Hq. 

Proof. It follows immediately from the spectral theorem and the fact that 
H{ is a nonnegative operator that, if ^ ) denotes a smoothed version 

of the characteristic function of the interval indicated, then 
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where 

hp := -p^ + Eix). 

Our aim is now to use the functional calculus for pseudodifferential op- 
erators developed in |HeRoj (see also |DiSj| , chapter 8) to show that, if 
^ e C^(]R), ^ = 1 on a neighborhood of [0, + E^], then 

CiK)xtlE^,.Jh,) = OienciJ^), (109) 

where hp := 

Once we have shown this, we will have 

but, applying lemma [HI we can write 

^Chp)xiho)^ = lim e(^)^n, 

n—>oo 

where the \E'„ are finite linear combinations of vectors of the form (11061) . with 

Lf = X(-oo,A)(/ip)v, < A < c^, SO that 

e(Mx(M^ = lime(M^'n= limxSol2c,+i?^)(Me(M^n 

= x^S,2c^+E^)^K)aK)x{h,)^ 

^ x{ho) = x[tc,+E^)Cho)x{ho) + 0{enc(^). 

We proceed now to prove (11091) . We first recall some facts we need from 
|HeRoj and |DiSj| . Given a function x ^ C*^(M) and a pseudodifferential 
operator P with symbol in a suitable symbol class (for our aims it is enough to 
say that this holds for both hp and hp), then also x{P) is a pseudodifferential 
operator, with symbol 

oo 

X{P) = Opf (a), ar^^e^aj, 

j=0 

2i-i ^ 

k=l 

where Op^ denotes the Weyl quantization, po is the principal symbol of P 
and the coefficients dj^k depend on the higher order terms in the expansion 
of the symbol of P (their precise form is given in |HeRo] ) . 
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We remark that the previous expressions are local in pq, and that for hp 
and hp the symbol is just the principal symbol, and is given by 

^o{x,p) := + E{x), t)o{x,p) := 

If we multiply two pseudodifferential operators, the symbol of the product 
is given by the twisted product of the symbols of the two operators involved: 

Opf (ai) ■ Opf (aa) = Opf (aitj.as). 



le 



(Vp ■ V^r - V| ■ Vg) ai(g,p)a2(x,0 



:iio) 



Applying these formulas to calculate the product i'^{h^)x^^lE^ c )(^p) ^^"^ 
using the locality in the principal symbol, we get that all the terms in the 
expansion of the product vanish, i.e., equation fllOQp . □ □ 



Corollary 4. Given a function x ^ C'o°(I^) a a > 0, we have 



\k\ 

where ^ G C, 
where 



oo / 
I 



:) and 



\k\ 



QMXiho) + Oois' 



2c^ + £"00 + niin{c^. A} 
sup{|A;| : k G supp x\ > 
sup , 



and we can choose sup{|A;| : k G supp ^} arbitrarily close to cg. 
An analogous statement holds for the creation operator. 

Proof. Applying twice lemma [HI the thesis will follow if we prove (11 111) re- 
placing /lo with ho. 

Applying lemma [HI we have 



X(-oo,Cx) 



{ho)QMX{hQ)'^ = QMX{ho)^ = lim 



where \E'„ is a linear combination of vectors of the form described in (11061) . 
with (f = X(-oo,A)(p^)v^- We have therefore 



idjV^{x, ■) 
1^1 



n— >oo 
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1^1 



^ n- 



Applying the operator to a vector of the form (11061) we get 
fid,v^{x,-)\ ^ idjv„{x,k) 



L2(R3,dfc) 



1=1 

■a^giY ■ ■ ■ a{gi)* ■ ■ ■ a^QMY^F 



with 



The functions gi have by hypothesis compact support in /c, with radius uni- 
formly bounded by c^. 

For the part depending on x, calculating the Fourier transform with the 
convolution theorem, we have 



:f( v{x)C^^^^^.9m)yp) = -j dk^J^^^^}^^g,(k).,MP + Kl 



where 

kj^ = (0, . . . , fc, 0, . . . 0) G M^^, 

with the k in the entry ji. Now, since Qj^, gi and (f have compact support, also 
the Fourier transform will have compact support in p, with radius bounded 
by + minjc;^, K} =: d^. This means that 



idjV„{x, ■) 
\k\ 



LP (g) 0(5-1)* ■ --aigMy^F 



idjV^{x, ■] 



= X{-oo,di){ho)ay — jip® aigiY ■ --aigMY^F, 

so, remarking that is a positive operator and that we can smooth X[o,d^) 
by a function with arbitrarily close support, we prove the thesis. □ □ 

Theorem 3. (Zero order approximation to the time evolution) The following 
two estimates hold: 

||(e-i^^I..| -e-"^oi)QMX{mres)\\ci.^^) = 0{VM+l\t\ey/\og{a{e)-^)), (112) 



||gM(e-^-., _e-'^o,)^(^.^^j||^^^,^ = 0{^/M+^\t\ey/\ogiais)-^)), (113) 
for every x G C^(M) such that XX = X- 
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Corollary 5. The subspaces associated to the Qms are almost invariant with 
respect to the dynamics generated by -ff^res ^ '^■^■> 

||[e-*^5../^,Q,,];^(if^^^J||^(^.^ = 0{VMTT\t\e^\og{a{ey^)). (114) 

Remark 6. (Adiabatic invariance of M -photons dressed particles subspaces). 

Using the unitary ^ we can translate all the previous results from the 
representation space to the original dynamics. This means that if we define 
the perturbed dressed projectors 

Pli := ^*QAf'^, (115) 

which satisfy by construction 

\\Pm - ^o'WciM') = 0{e^logia{e)-^)), 

we get that 

\\[e~-^'''''"/^,PlMH'n\\ci.^^) = 0{VMTl\t\eV\ogiaie))-^), 

and 

||^^_i,H-/e^^M] ~(^.,.)||^^^^^ = 0{{l + \t\)e^\ogiais)-^)), Vx G C^iR) . 

For the original dynamics we have then 

[e-'''''/',P!,]xm = [(e-"^^/^ -e-'*^^"/^),PX,]x(i^^) + 

+ [e-''^'^^/^Pl,]{xm - xiH^n) + [e-"^^"/^P|,]x(i^^''^) . 

The first term can be bounded by 0{a{ey^'^e^^) = 0{e'^) using proposition 
[Hand equation ([HHD- The second one by C((T(e)^/^) using lemma El and the 
third one has been just estimated above. 

Putting together these facts, equation ( |2T]) is proved. 

For the particular case M = 0, this result, together with the expression 
of the zero order Hamiltonian Hq, was already shown in [Te^j, assuming an 
infrared regularized interaction and a relativistic dispersion relation for the 
particles, which automatically implies that they have a bounded maximal 
velocity, and avoids therefore the introduction of cutoff functions. 

Proof. First of all we remark that, employing lemma [3, we can replace 
x(-ffdres) with x(^o) up to tcrms of order 0{e)c(._^), which are smaller than 
the error we want to prove. 
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Since H^^^^ and Hq are both self-adjoint on Mq, we can apply the Duhamel 
formula, and obtain 

= [ ^ise^(^-*)^-==/^(i/^,,3-/zo)e-^^^»/^gMx(/^o) 
^ Jo 

Jo 

= -i [ cise^(^-*)^S.os/-/,,_^Q^,^(/,o)e--'^o/. + 0{eVMTlV^og{a-^)) , 
Jo 

which implies 

< \t\ ■ \\hi,^QMX{ho)\\ciM^) + 0{e\t\VMTlV'^og{cr-^)). 

We proceed now as in the proof of lemma [71 From equation (!95l) it follows 
(omitting the terms of the form "p ■ $" , which can be treated in the same 
way) 



- ^(1 - x{V:H''''V,))<!>\\V,v,) ■px{V:H'''^V,] 



QuXiho) 



= \^iV,v„)-pQMX{h) - \xiy:H'^''V,)^'{iV.v,)-pQMX{ho) 

- i(l - xiy:H''''V„))^'{iV,v„)-pQMX{ho) 

- i(l - x{V:H'^''V,))^'{:iV^v,)-p\x{V:H'^''V„) - x(/^o)]Qmx(/^o) 
= - x{hom\^,v^) ■ p[x{V:H''''V^) - x{ho)]QMX{ho) 

+ 1[x{V:H''''V^) - x(/^o)]$'(iV.t;.) ■ p[x{V:H''''V^) - x{ho)]QMX{ho) 

Applying lemma [71 we have immediately that the last term is of order 
We examine therefore more closely the second one. 

Using equation (11021) and the corollary to lemma [HI we get that 

[x{V:H'^'^V^) - xiho)]QMXiho) = 
= (Qm+1 + ClM-i)aho)[x{V:H''''V,) - x{ho)] ■ Quxiho) + 0(^')z:(^) , 
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where the function ^ has a support shghtly larger than that of the function 
X- Applying again the corollary, we get that, if we choose the support of the 
function x sufficiently smaller than the support of the function x, then 

hi,^QMX{K) = 0{e^), (116) 

so we get the first estimate. 

For the second one, we apply again the Duhamel formula, but inverting 
the position of the two unitaries: 



drcs/ 



e 



- Jo 

= -i [ rfse'(^-*)'^°/^QM/ii,^x(^dreJe-^^^S-/^ + 0(£VMTTVlog((T-i)) . 
Jo 

It follows from the proof of lemma [7] (in particular equation (11041) and what 
follows) that 

QMhl,^X{mres) = 0(£v/MTTVl0g(0--l))^(^), 

SO also the second estimate is proved. □ □ 

Lemma 10. The truncated dressed Hamiltonian 

H!il,^:=ho + eh,^ + e'h,^ (117) 
is self-adjoint on J^q for e small enough. 

Proof. By construction, the coefficients hi^^ belong to C{J^, J^), and define 
symmetric operators on JifQ. Moreover Hq is self-adjoint on Jifo, therefore 

ii(^dl,x - homu < csiwhon.-^^ + w^iu), g .^o. 

By a symmetric version of the Kato theorem ( |ReSi2| , theorem X.13) the 
claim follows. □ □ 

Theorem 4. (First order approximation to the time evolution) Given a 
function x G 



dres7 



Jo 



+ 0{e'/'\t\)ci.^^){l - 6mo) + 0{e\\t\ + \t\') ^/\og{a-^)) ^^^^ , (118) 
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where 6 mo = 1 when M = 0, otherwise, [Hj^ ,Qm] = WM, 
^ 1 

■= ll^eP" + E{x)+H, (119) 
1=1 ' 



^(2) 



Jl=l 



2 

wt/i mf := 1/(1 + yC/), K := k/\k\ and 



er.= - [ dx dy '-P^ . (120) 



The off-diagonal Hamiltonian is defined by 

h2,OB ■■= $a • VE(X) . (121) 

Proof. We split the proof into three parts. In the first one, we show that 
equation (11181) is true with a diagonal Hamiltonian H^^ given by 

Af 2 

1=1 j=l 

+ («:-pJ(«:-p,.Je'^-(-^i--'i)], 
and an off-diagonal one /i2,OD defined by 

^2,OD := $a ■ VE{x) - - ^{dj^o^ipipj + pipjdj^o 

^ 3N 

In the second part we prove that if one neglects the term 

— ^ a{idjV„ya{idjVa-) 

A — 1 



(123) 



_2 37V 
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in one gets an error of order in the time evolution. Note that 

this term is exactly zero if the initial state for the field is the Fock vacuum. 
In the third part, we prove analogously that we can replace /i2,oD with /i2,oD- 

More specifically, the terms which we neglect in Hj^ and /?.2.od give rise 
to higher order contributions to the time evolution, although their norm in 
C{J^o,J^) is not small. This is caused by the fact that they are strongly 
oscillating in \k\, so that their behavior is determined by the value of the 
density of states in a neighborhood of = 0. For all these terms, the density 
however vanishes for k = 0, uniformly in a, and this implies that they are 
of lower order with respect to the leading piece $o- ■ V-E(x), whose density 
instead diverges logarithmically in a (we elaborate on this last observation 
in a corollary to this theorem). 

We start showing that we can, up to the desired error, replace -ff^res ^y 

^dies,x- lemma fTOl -f^dres,x self-adjoint on J^q like -ff^res; therefore we 
can apply the Duhamel formula and use theorem [2] to get 

e-"^S.esA _ e-<L,x/- = -1 rrfse'(^-*)^-=/^(ffd',,,-i/ilje-^^^-- 

^ Jo 

Moreover, using lemma [71 we can replace QMX(-f^dres) by Q MX{ho)x{Hl^^^) . 

Since the diagonal Hamiltonian H-^ is also self-adjoint on for e suffi- 
ciently small (the proof can be given along the same lines of lemma [TU|l . we 
apply again the Duhamel formula, 

^ Jo 

= -i rfse^(^-*)^--x/^/ii,^e-'^</^QMX(/io)^(i^^ 

^0 

-ie f ds e^(-*)^^l,x/^(/,2,^ - h^,^)e-'^<''QMx{ho)mirJ 
Jo 

where 

3Af 



^dres ) 



/i2,D ■.=^^a{\djV^)*a(\djV^) 

7 . „■ . 1 J I I 

+ ('^■Pj('^-p,je''-^^^i-^''i)]. 



CiT^ii) 
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To analyze the first term, we remark that, proceeding as in lemma [3, one 
can prove that 

x{ho) - x{hS^) = 0(eVlog(a-i))£(^), 
so 

therefore, with equation (11161) . 

^0 

= _i rdse'(^-*)^-=.x/-/ii,^Q,,^(/^o)e-'^</^^(i/^,,J 

+ 0{e'\t\^logia~^))ci.^)=0{e'\t\^\ogia-^))ci^). 

Concerning the second one, applying once again the Duhamel formula, we 
have 

-ie / ds e'(-*)^il.x/-(/i2,^ - h2,^)e-'<''/^QMX{homHLs) = 
Jo 

= -is [ dse'^'-'^''^''^^J'{h2,,,-h,T,)e-'''''^ 
Jo 

+ 0(£2|t|Vlog(^-^))£(^), 

so we have to look at e'*-'*~*^'^drcs,x/^(/i2 — /i2,d)Qa/x(^o)- Following a pro- 
cedure already employed several times, we first observe that, in the expres- 
sion for /z2,x, equation ( |971) . we can replace, making an error of order 0{e), 
xiy:H^^-V„) with x{ho). 

Using the corollary to lemma [HI we can then eliminate from /i2,x'5A/x(^o) 
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all the terms containing (1 — x(^o))- What remains is then 

(/l2,x - h2,T>)QMX{ho) = 

2 r 1 



4 

i=-2 ^ j=l 



2 

1 



2 

i=-2 



+ O(£Vlog(0r-l))£(,^). 



(2) 

Applying theorem [HI we can now replace e'*-*~*^'^drcs,x/^ with Q^i^-'^)^^/^ making 
an error of order 0{\t\e^J\og{a-^)). 

Concerning the last term, we get, expanding the commutator and apply- 
ing again the corollary to lemma [9], 

1 

= Yl QM+.e(/^o) [^\iV,v„) -p, $J • p]QMX{ho) 

i=-l 
1 

= Yl QM+^^{ho) [^'{id^v,), K,i]PiPjQMX{ho) + 0(£Vlog(a-i))£(^) . 

i=-l 

Using equation fl96p . we can check that this cancels exactly the diagonal 
terms in h2^x ~ ^2,d, so that 



= e'^'-'^''^'/%oDQMXiho) + 0{e\t\ v^M^) . 

~ (2) (2) 

We proceed now to show that we can replace with , up to an 
error of order 0{e'^^'^\t\)jr{jf)- 

Applying repeatedly Duhamel formula, and putting \& := Qmx{HIj.^J'^q, 
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we get 

3N „t 



4 Jo 

+ 0{e'\t\')ci^y 



To streamline the presentation, we assume that M = 1, the calculations for 
M > 1 are basically the same, but more cumbersome. 
The time integral gives 



^g^is{\kx\-\k\)/£^ishp/e^iXj\ki-k)^-ishp/e-^(^^^ = 

l^il^/^^ |A;|V2[l + i(|A:i|-|A:|).-i]^^ + ^^'^^' ^ 

^ggis(lfcil-|fcl)/egis/ip/egixj-(fci-fc)g-is/ip/e^^^^ 

(124) 

Integrating by parts we get 

i(|A;i| - \k\)e-^ [ rfse^'(l'=^l-l'^l)/V^p/V^^-('=-'^i)e-^''^p/"^(x,A;) = 
Jo 

— QH\ki\-\k\)/e^ithp/e^iXj-{k-ki)^-ithp/e.qj _ giaJj -(fc-fci) ^_|_ 



- [ c/se''(l'^il-l^l)/V'p[/ip,e'^'^-('=-'=^)]e-^''*''^, 
^ Jo 



where the commutator is of order 0{6) when applied to functions of bounded 
kinetic energy, so that the right-hand side is uniformly bounded in e. 

We have now to put this expression back in (11241) and estimate the single 
terms. We show how to do this for the first one, the others being entirely 
analogous. We ignore the unitary on the left, which does not change the 
norm, so we have to consider 

0a{kl)_^2 f ■ l^)^o{k) 



\kmi+m\-\k\)e-^] 

QH\ki\-\k\)/£^ishp/e^ixj-(ki-k)^-ishp/e.qj^^^ 
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for one fixed j. Using twice the Cauchy-Schwarz inequality we get (we put 
for brevity f{k) := ejip{k)\k\-^/^) 



|2 



^ [ dx I dki\---\^< \t\ I dx I dkidk\f{ki)\^- 

Jr^N J^3 J J 1 



+ m-\k\rs- 



dk / ds 



gis/ipgixj-(fc— fci)g— is?ip 



<Ce'mni, / dk. 



A/e 



A/e 



dk 



kik 



^0{e\t\'\mlf)- 



/o 1 + (fci - k)^ 
We examine now separately the last two terms in ^2 od- For the first we 



get 



ie f ds e^(^-*)'^°/^a^.$o,iMe""''°/'QMX(i^dres)^o 
^0 



V2 



— — le 



—itho/s 



dse'''"'/' 



didjv{x, k) 



+ a 



didjv{x, k) 



\k\ 



PjPie 



V2J0 

The part with the annihilation operator gives 



|/t|V2 J 



|A;|V2 ; 



V2 



I^32'^l2^ 1/2 



-itho/e 



V2 



.g-i.ftpAg-i.|fe|/eg-i.E;r=-i'l'=Ml/^^(a;. k,ki,..., ku-i) 

V2^ L \k\'/' Jo ^ ^''^ 



V2 U \k\'/' 
■e-''\^\l'-^{x-kM,---.kM-i 



g-it^oA^^ / dk 



\k\y\l-i\k\e' 



— (l-i|A;|£-') / dse''^^'' 
) Jo 



2 -ishp/e -is\k\/£ 



-^{x;k,ki,...,kM-i) ■ 



(125) 
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Integrating by parts we have 



Jo 

Jo 

= e^ih^/^e--^k-x,^ (fi; . p^.j2e-i*/^p/^e-'*l'=l/"^(x; k,ki,..., ku-i) 
- e-'^^-^w {k ■ k,ki,..., ku-i) 

_ 1 /" dse''''^''[h^,e-'^-''^^{K ■ p^.j2]e-'^'^p/"e-^^l^'/"^(a;; k,ki,..., ku-i) 
^ Jo 

= e^*'^p/"e-''=-"^i(/€ ■ p^.j2e-'*'^p/"e-'*l'=l/"^(a;; k,ki,..., ku-i) 



t 



-1 / dse''^^'' 
'o 

-ik-x 



-ik-Xi 



[K-p 



-''''-''e-''\^\''^{x-kM.---Mi-i) 



^tK/e^-ik-x,^ ■ p^.j2e-'*'^p/"e-'*l'=l/"^(x; k,ki,..., kM-i) 



- e-"=-^« {k ■ Pj^^ix; k,ku..., kM-i) 



+ i / dse''^-^!' 



e ''"■^^K . p.^ - £\k\e 



[K-p 



■ji' 



We have to put now this expression back into (11251) and estimate the result. 
We show how to proceed for the most singular term, i. e. the one containing 

i ! dse''^^'%-'^-''^^{K ■ p^.j3e-'^'^p/"e-^^'^l/"|fc|^(x; k,k^,..., ^m-i), 
Jo 

the others can be treated in the same way. Putting this term back in (11251) . 
and ignoring the unitary e"'*^"/"^, which does not change the norm, we have 
to estimate in the end 



i^^M f dk 

V2 h.\k\y\l-me ,Jo 



- / rfse'^'^p/"e-''=-"^i(/€-p^.jV 



ishp/e _ 



■e 



--l^l/^|fc|vl/(x;A;,fci,...,A:M-i) 
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Applying Cauchy-Schwarz inequality we get 



|A;|(1 + |A;|%-2) 

■'^{x;k, A;i,...,/cm-i)'^ 
The first integral gives 

dk,^JM$L_^ ^ c f\k — ^ = cs' ^ 



= 0(£2l0g(l/£)), 

so, calculating the norm we get 

I- ■ -115^ < C£^log(l/£)|t|'M [ ds ||e'^'^p/^e-^'=-^^i- 

Jo 

■{k . p^J'e-'^^-/^e-'^\'\/^\k\^l'{x; k,k,,..., ^m-i) ^ 

= CeHog{l/e)\t\''M f ds\\{K ■ p^Jh-''''-/'\k\^{x;k,ki, . . . ,kM-i 

Jo 



which shows that the norm of this term is of order 0{e'^\t\ y^log{l/e)) in 
C{J^). The part with the creation operator gives 



— J- 



V2(M + 1) |A;JV2 
•^'(a;;/ci,...,A;^,...,A;M+i) 



.y2{MTT)\k,mi-i\k,\e-^r ' Vo 
(«;^ • p,.J^e-^^'^-/^e-'^l'=''i/^^(a;; fci, . . . , . . . , ^m+i) • 
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We can now integrate by parts using equation fll26p . As for the annihilation 
part, we examine just one term, 



- [ c/se'''^p/V^^-^^i(/t^-p^J^e-^''^p/"- 
) Jo 



v/2(M + l)(l-i|A; 



The norm squared is given by 



2(M + 1) 



dx / dki . . . dk 



M+l- 



\Qh{kf,)\^\k^ 



P,fe-'^'^^/'^{x-k,,...,k„...,kM+i: 



2 



e 



2(M+ 1) 



dx I dki . . . dkM+i 



\Qn{^k^)\^\k^, 



dg Q-'^s\k^,\^ishp/e^iekf,-Xj 



P,fe''^'^^/'^{x-k,,...,k„...,kM+i: 



< 



e'\t\ 
2(M+ 1) 



dki ■ ■ ■ dkpf+i 



e'\t\ 
2(M+ 1) 



Pi,fe ''^^/'^{x; fci, . . . , fc^, . . . , kM+i 
\gj,{ekf,)\'^\kf,\ 



dki ■ ■ ■ dkjif+i 



^^^^ ■ Pj,fe-'''''^/'^{x; fci, . . . , A;^, . . . , kM+i] 



ds ||e ^^\^^\ Q^^^p/^ (^^^^^^'"^31 . 

2 



e'\t\ 



2(M+ 1) 



dki ■ ■ ■ dk^j^i 



■"^{x^ki, . . . ,kf„ . . . ,kM+i] 



\gj,{ek^)\^\k^ 



,dx) 



ds\\{K^-p^ye-''^^/'- 



< 



2(M+ 1) 



00 /£ 



d\ki, 



■'^{x;ki,...,k^,...,kM+i 

Ce'\t\gl 
4(M+ 1) 

■QAfX(-f^drcs)*o||^, 



Ml 
i2 



dn^ / ds \\{k^ ■ p j e 



3 —ishp/e 



- log(£-2 + 1)] j dQ^j^ ds\\{t,,-p^fe-''^^/' 



so we get that this term gives a contribution of order in the norm of 
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We separate also in the second term of /i2,oD the annihilation and the 
creation part. For the annihilation part we get 



le 



3N 



3=1 



^ _^^M(M-l)e-'*''«/^ J] / di / dUdn 



dg g-i«(l?l + KI)/£gisV/£g-i(€+C)-2;jig-is/ip/£^ _ 



We proceed now in the same way as we did for the first term. Integrating by 
parts we get 

-i(iei + ici) 



— g-i*(ICI + ICI)/egit^p/£g-i{?+C)-a:jig-it/ip/ev]> — g-i(5+C)-a:ji v|/ 



As in the previous case, we examine just one term, 



le 



^M{M-l) dU dCi,,Q,-— 



ie|V2|c|V2[i-i(|ei + ici)£-i] 

• / ds e~'**^l^l+l^l)/V'''^p/^e'"^^^"^^^'^^i + C) ■ P ■ e'"'*'^p'^^ 
Jo 



the others can be treated in the same way. Using Cauchy-Schwarz inequality 
we get that, independently of a, 



l^<Ce'\t\M{M-l) d\i\j^ d\C\^ 



l eilCKieP+ICP) 



ds 



l^ + CI 



l^oll^ 



4|j.|2\ 



The creation part can be estimated in a way entirely analogous to the 
one already employed for the first term of ^2,00- n □ 



57 



Corollary 6. The radiated piece (i. e. the piece of the wave function which 
makes a transition between the almost invariant suhspaces) for a system start- 
ing in the Fock vacuum is given by 

^rad(t) := (1 - go)e-''^5.os^(a;)(]p = 

__t^^_ii;,„/,^^j_i_;^_y dse^mi-O^Y {x''^{s-x,p))ij{x)+ (127) 
where 

\\R{t,e)U < Ce'logie-'m + + \\\x\^\U + \\\mU) , 

and x'j is the solution to the classical equations of motion 
x''j{s;x,p) = -V^^E{x%s;x,p)), 

(128) 

x''j{0;x,p) = Xj, x'j{0;x,p) = Pj, j = l,...,n. 

We get the leading order of the radiated piece corresponding to the original 
Hamiltonian H'^ , for a system starting in the approximate dressed vacuum 
^a{e){.x), applying to this wave function the dressing operator Va-(s){x). 

Proof. Applying equation flll8p for the case M = we get at the leading 
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order 

-t 



Jo 

v2 |fc| / Jo 



+ 0{e^\t\)ciM')U\\LHR^-), 

where we have used Egorov's theorem to approximate e'**''p/'^Va;^-E(x)e~''*'*p/^ 
(see, e. g., |Roj ). To end the proof we have to show that the norm of the 
first term is small. The procedure to employ is identical to the one applied 
several times in the proof of theorem HI First integrate by parts with respect 
to s and then estimate the resulting terms. For the sake of completeness we 
show how to estimate one of these terms 



v^|A;|3/2(l + i|A;|£ 
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because the others have an analogous structure. Its norm satisfies 
II l|2 



< 



dk 



mi + m Jo 



<Ce^\t\ / d\k\ 



A/e 



\k\ 



ds 



ds 



L2(R3JV) 



L2(R3iV) 



Ce^\t\ sup \W^^E{x)\\og 



<Ce'\tm + \t?)^og 



1 + 

" -IP, 



jajjie 



-is/ip/e 



L2{R3JV) 
2\ 



,l + a(£)2 

where in the last inequahty we have apphed theorem 2.1 from jRaSi J. □ □ 
Remark 7. The norm squared of the leading part of the radiated piece is 



I dx I dk 

371 



> 



e 



-2 



> 



e 



127r4 



dx I dk — 

k 



dx 



Opf dse-''^''Y^^x]{s-x,p)^i){x) 
Opf £ ds e''^'/^ ^' P)) ^(^) 

Opf Q*dse-'^'=f;i,^(s;x,p))^(, 



> 



C3 



> 



'a(e)e-i 

The symbol which appears in the Weyl quantization 

AT 



C3 



/ dse "'^y~^a;^(g;x,p) 
•^0 ,=1 



is independent of e and for /c = is a non null function, 

N N 

Y[x%t;x,p) - x'j{0,x,p)] = Y^[x''-{t;x,p) - p_ 



therefore the corresponding operator will be also non null. Therefore, we 
expect for a generic state ip that 



inf 

0<fc<A 



Opf £ ds e-'^' J2 ^"M^ ^' P)) ^(^) 



>0. (129) 



L2(lR3iV)0C3 



60 



Then the norm of the radiated piece will be bigger than 

TV 



inf 

0<fc<A 



I A I dk — 

K^'^T^ Ja(e)e~^ k 



1/2 



which gives a lower bound on the transition almost of the same order of the 
upper bound. 

Remark 8. The radiated energy, defined in equation flHTj) . can be written at 
the leading order as 

^rad(t) = (e-'*^d.eeA^,(a;)^]p, (1 ® i/f)e-"^d.e.A^,(a;)^]p), 

where ipd is defined in ( l30l) . Using the expression for the radiated piece we 
get 



N 



ds / rfs'e*(^-^')l''l/^ 



-y 

■(k ■ Opf (i^(s'))V^d, ■ Opf (^j'(s))7/;d)L2(R3iV) 



127r2 



> ejC,- ds ds' - 

h 'Jo Jo i( 



s — s 



,i{s-s')A/e _ 



1) 



■(^d,Opf (i-(s') ■ i^(s))^/^d)L2(K3iV) 



where we have inserted the explicit expression of the form factor given in ([3]) 
and used the product formula (11 101) at leading order. The radiated power is 
then 



Prad(t) - |i?rad(t) - |^ 



N 



j2 r^.^i^t^^-^)^/^] 



t - s 

which converges formally to the expression given in flH^ when e ^ 0~ 
Corollary 7. Let 

uj{t) := e-"^d.e./^cuoe"^d.e./^, 
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where ujq G J^i{Q MX{H^res)'^) > Banach space of trace class operators on 
Q MX{H^res)'^ ! ^'^^ j'arfza/ trace over the field states 

ujpit) := tr,j^ uj(t), 

then 



+ 0(£2|t|v/log(a(£)-i))^,,(i2(ij3^))+0(£2|t|Vlog(a(£)-i)) .i?fl(L2(]R3JV-)), 
Z=l ' 

and J^i(L^(]R^^)) denotes the space of trace class operators on L^(M^^). 
Proof. The proof follows from the following three facts: 

1. the term of order e in equation flllSp is off-diagonal with respect to the 
Qm] 

2. the diagonal Hamiltonian Hf^ , defined in f lllQp . is equal to -f^D.p ® 1 + 
1 (S> -fff, so we have that 

tr^(e-'<'/^a;oe'<'/^) = e-'<p/^tr^(^o)e"<p/^ ; 

3. the following well know inequality, which holds for any Hilbert space 

□ 

Theorem 5. Given a macroscopic observable for the particles, Op^(a), 
where a is a smooth function bounded with all its derivatives, and a den- 
sity matrix u G J^i{Pljx{H'^)J^) whose time evolution is defined by 

uj{t) := e-'*^'/^we'*^'/^ (130) 
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then 

tr^(^(Opf (a) ® l^)oo{t)^ = tr^2(M3.)(^Opf (a)e-'<itr^cD(0)e"<pJ + 

- 5mo) + 0{e\\t\ + |tn log(a(£)-i)), (131) 

where 

m := . (132) 

Proof. First of all we observe that, using proposition [Hand lemma [3l we have 



tr^-(^(Opf (a) ® l^)a;(t)^ = tr^(^(Opf (a) ® 
■ e'*^^'"'^'/.^ + 0(^(5)^/') + 0{a{e)e-') , 



'UJ. 



where uJa{£) ^ "^iiPMX{H'^''^^'^^)J^) ■ By the definition of the dressed Hamil- 
tonian and the cyclicity of the trace we have then at the leading order 

tr^(^(Opf (a) ® l^)u;{t)^ = ^(^(e)^^) ^ oia{e)e-') + 

+ ti^' (Opf (a) ® l,^)^*e-'*^drcB/^<^^^(^)^*e'*^d.csA^ . 

The transformed observable, using the definition of and lemma [TJ is given 

by 

^(Opf (a) ® 1^)^* = V:{OvY{a) ® 1^)V^ + iexiho)^i ■ pV; ■ 

■(Opf (a) ® + i^l - xiho))K-pxiho)V:{OpYia) ® 

-i£V;(Opr(a) ® U)VMho)p ■ K 

-i6V:{0pY{a) ® 1^)K(1 - x(/^o))p- $^x(/^o) + 0(£'log(a-i)) . 

The operator can be considered as the Weyl quantization of the operator 
valued symbol 

(a;, p) ^e-'*^'"'^ G C{^) . 

It is not in the standard symbol classes because the derivative is an un- 
bounded operator, but if we multiply it by Qm we get a smooth bounded 
symbol. In the calculation of the trace is always multiplied by Qm or 
Qm±i, therefore we can proceed as if it were in a standard symbol class. 

A simple application of the product rule for pseudodifferential operators 
(equation (A. 11) [TesQ gives then 

V:{OpYia) ® = Opf (a) 1,^ + V;*[(Opf (a) ® 1^), K] = 

= Opf (a) ® 1^ - ieV:OpY{{a, V^}) + 0{e')ci.^^) , 

63 



where {■, ■} denotes the Poisson bracket. We get therefore 

{a,K} = iVpa)-i^iV.v^)V^ 

and 

r;(Opf (a)® 1^)V; = Opf (a)®l^ + e<l>(iV.t;.)(Opf(Vpa)® 1^) 

+ 0{eHog{a-')). 

Using this expression we have 

'^(Opf (a) ® 1^)^* = Opf (a) ® 1^ + £$(iV.t;^)(Opf (Vpa) ® U) + 

+ iex{ho)<^i ■ p(Opf (a) ® 1,^) + i£(l - x(/^o))$;^ ■ PX{ho) (Opfla) ® 1^) 
-i£(Opf(a)®l^)x(Mp-$^ 

- i£(Opf (a) ® 1^) (1 - x{ho))p ■ <^ix{ho) + 0{e' log{a-')) . 

All the terms of order e in the previous expression are off-diagonal with 
respect to the Qjv/S, and the same holds for the term of order e in (11181) . 
Therefore, they all vanish when we calculate the trace. Using point 2 and 3 
of last corollary we get then (11311) with 

Using again lemma [3] and the fact that the terms of order e in the expansion 
of ^ are off-diagonal we get in the end also (11321) . □ □ 
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